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Abstract 

We propose a two Higgs doublet Type III seesaw model with /i-r flavor symmetry. We 
add an additional SU(2) Higgs doublet and three SU(2) fermion triplets in our model. 
The presence of two Higgs doublets allows for natural explanation of small neutrino 
masses with triplet fermions in the 100 GeV mass range, without flne tuning of the 
Yukawa couplings to extremely small values. The triplet fermions couple to the gauge 
bosons and can be thus produced at the LHC. We study in detail the effective cross- 
sections for the production and subsequent decays of these heavy exotic fermions. We 
show for the first time that the ^.-t flavor symmetry in the low energy neutrino mass 
matrix results in mixing matrices for the neutral and charged heavy fermions that 
are not unity and which carry the flavor symmetry pattern. This flavor structure can 
be observed in the decays of the heavy fermions at LHC. The large Yukawa couplings 
in our model result in the decay of the heavy fermions into lighter leptons and Higgs 
with a decay rate which is about 10^^ times larger than what is expected for the 
one Higgs Type HI seesaw model with 100 GeV triplet fermions. The smallness of 
neutrino masses constrains the neutral Higgs mixing angle sin a in our model in such 
a way that the heavy fermions decay into the lighter neutral CP even Higgs h^ , CP 
odd Higgs AP and the charged Higgs , but almost never to the heavier neutral CP 
even Higgs . The small value for sin a also results in a very long lifetime for /i". 
This displaced decay vertex should be visible at LHC. We provide an exhaustive list 
of collider signature channels for our model and identify those that have very large 
effective cross-sections at LHC and almost no standard model background. 
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1 Introduction 



Understanding of the flavor structure of tlie fermions has emerged as one of the most 
formidable problems in particle physics. While all fermions are expected to attain masses 
in the standard model through their Yukawa couplings with the standard Higgs doublet, it 
is not clear why the mass of the electron should be six order of magnitude smaller than that 
of the top quark. The extremely tiny neutrino masses pose a further challenge and demand 
an explanation. The "seesaw" mechanism [1] has been the most widely accepted method 
of explaining the smallness of the neutrino mass compared to that of the charge leptons. 
The seesaw mechanism was so named because the lightness of the standard neutrino is 
explained due to the heaviness of an additional particle beyond the standard model of 
particle physics. The new mass scale could be associated with a GUT scale, or in general 
with any intermediate mass scale. Being much heavier than the rest of the standard model 
particles, this additional field can be integrated out, giving a dimension five Majorana 
mass term for the neutrinos [2]. This mass term is inversely proportional to the mass of 
the heavy particle, and hence neutrinos become naturally light. There are three variants 
of the seesaw mechanism. These come from the fact that one can obtain the dimension 
five effective operator by integrating out either SU(2)iX U(l)y (i) singlet fermions [1], or 
(ii) triplet scalars [3], or triplet fermions [4] (also [5]). The three variants are commonly 
known as Type I, Type II and Type III seesaw mechanism, respectively. While Type I and 
Type II scenarios have been extensively explored in the literature for a long time, focus 
has only recently shifted to the Type III seesaw mechanism, and a plethora of papers have 
appeared of late. While the possibility of gauged U(l) symmetry with fermion triplets 
was studied in [6], authors of [7] studied for the first time predictions for leptogenesis 
within the framework of Type III seesaw. A hybrid Type I+III seesaw framework is shown 
to result within a SU(5) GUT model in [8.9]. and within a left-right symmetric model 
with spontaneously broken parity in [10]. The effect of the additional fermions on the 
Higgs mass bounds was studied through renormalization group equations in [11], while the 
renormalization group evolution of the neutrino mass matrix within the Type III seesaw 
framework was performed in [12]. In [13] the authors work with just one extra heavy 
fermion triplet and generate the addition light neutrino masses at the loop level. The 
phenomenology of the Type III seesaw in lepton flavor violating processes was studied in 
great depths in [14,15] and also recently in [16]. 

The most crucial feature concerning the Type III seesaw is the following. Since the 
additional heavy fermions belong to the adjoint representation of SU(2), they have gauge 
interactions. This makes it easier to produce them in collider experiments. With the LHC 
all set to take data, it is pertinent to check the viability of testing the seesaw models at 
colliders. The implications of the Type III seesaw at LHG was first studied in [17] and [18] 
in the context of a SU(5) GUT model. In the SU(5) model it is possible to naturally have 
the adjoint fermions in the 100 GeV to 1 TeV mass range, throwing up the possibility of 
observing them at LHC. The authors of these papers identified the dilepton channel with 
4 jets as the signature of the triplet fermions. Subsequently, a lot of work has followed on 
testing Type HI seesaw at LHC [19-21]. 
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In the usual Type III (and also Type I) version of the seesaw models with one Higgs 
doublet, the neutrino mass is given by 

m. = -v'Yi^Yj: (1) 

where, v is the Higgs Vacuum Expectation Value (VEV), M is the mass (matrix) of the 
adjoint fermions and Y-^ is the Yukawa coupling (matrix) of these fermions with the stan- 
dard model lepton doublets and Higgs. To predict neutrino masses ~ 0.1 eV without fine 
tuning the Yukawas, one requires that M ~ 10^^ GeV. On the other hand, an essential 
requisite of producing the heavy fermion triplet signatures at the LHC, is that they should 
not be heavier than a few hundred GeV. One can immediately see that if M ~ 300 GeV, 
then nil, ~ 0.1 eV demands that the Yukawa coupling Y-^ ~ 10~^. This in a way tenta- 
mounts to fine tuning of the Yukawas, and smothers out the very motivation for the seesaw 
mechanism - which was to explain the smallness of the neutrino mass without unnaturally 
reducing the Yukawa couplings. 

In this paper, we propose a seesaw model with 300-800 GeV mass range triplet fermions, 
without any drastic reduction of the Yukawa couplings. We do that by introducing an 
additional Higgs doublet in our model. We impose a Z2 symmetry which ensures that this 
extra Higgs doublet couples to only the exotic triplet fermions, while the standard Higgs 
couples to all other standard model particles [22] . As a result the smallness of the neutrino 
masses can be explained from the the smallness of the VEV of the second Higgs doublet, 
while all standard model fermions get their masses from the VEV of the standard Higgs. 
Therefore, we use the presence of two different VEVs in our model to explain the smallness 
of the neutrino masses compared to all others, without resorting to unnatural suppression of 
the neutrino Yukawa couplings. We show that these large Yukawas result in extremely fast 
decay rates for the heavy fermions in our model and hence have observational consequences 
for the heavy fermion phenomenology at LHC. We show how this can be used to distinguish 
our two Higgs doublet Type HI seesaw model from the usual one Higgs doublet models|3 

The presence of two Higgs doublets in our model also enhances the richness of the 
phenomenology at LHC. We have in our model two neutral physical scalar and one neutral 
physical pseudoscalar and a pair of charged scalars. We will work out in detail our Higgs 
mass spectrum by imposing constraints coming from the neutrino masses. We will show 
that due to these constraints, our Higgs mixing angle is very small and the Higgs behave 
in a very peculiar way and have collider signatures which are very different from the usual 
two Higgs doublet models in the market [23-27]. We will study this crucial link between 
neutrino and Higgs physics in our model and its implications for LHC in detail. 

Another feature associated with neutrinos which has puzzled model builders is it unique 
mixing pattern. While all mixing angles are tiny in the quark sector, for the leptons we 
have observed two large and one small mixing angle. In its standard parametrization with 

"'^The largeness of the Yukawa couplings (along with the smallness of the heavy fermion masses), also 
brings in larger non-unitarity and larger lepton flavor violation in our model, compared to the earlier Type 
III seesaw models. However, these are still well below the sensitivity of the current and upcoming future 
experiments. 
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mixing angles Ou, O23 and ^13 and phases 5 (Dirac), a and (3 (Majorana), the neutrino 
mixing matrix is given as 

C12C13 S12C13 si3e~'^ ^ 

UpMNS = -C23 •S12 - 523 Si3 C12 e"^ C23 C12 - S23 Sl3 Sl2 523 C13 

\ S23 S12 - C23 Si3 C12 e'^ -S23 C12 - C23 Si3 S12 e'^ C23 Ci3 / 

In this parametrization, the mixing angle ^23 is observed to be very close to 7r/4, while 
^13 has so far been seen to be consistent with zero. This indicates that there should be 
some underlying symmetry which drives one mixing angle to be maximal and another to 
be zero. The most simple way of generating this is by imposing a /x-r exchange symmetry 
on the low energy neutrino mass matrix [28] . 

In this paper we will impose the fi-r symmetry on the Yukawa couplings and the heavy 
fermion mass matrices. This leads to fi-r symmetry in the light neutrino mass matrix and 
hence the correct predictions for the neutrino oscillation data. We discuss in detail the 
light as well as heavy neutrino mixing. We first provide general expressions for all mass 
eigenvalues and mixing matrices and then study the experimental consequences for our 
model. We show that due to the fi-r symmetry, the mixing matrices of the heavy fermions 
turn out to be non-trivial. In particular, they are also /i-r symmetric and hence much 
deviated from unity, even though we start with a real and diagonal Majorana mass matrix 
for the heavy triplets. This is a new result and we will show that this affects the fiavor 
structure of the heavy fermion decays at colliders, which can be used to test //-r symmetry 
in neutrinos at LHC We study in detail the collider phenomenology of this fi-r symmetric 
model with three heavy fermion SU(2) triplets and two Higgs SU(2) doublets and give 
predictions for LHC. 

The paper is organized as follows. In section 2, we present the lepton Yukawa part of 
the model within a general framework and give expressions for the masses and mixings 
of the charged and neutral components of both light as well as heavy leptons. In section 
3, we present our fi-r symmetric model and give specific forms for the mass and mixing 
parameters. We show that the mixing for heavy fermions is non-trivial and /x-r symmetric. 
In section 4, we study the cross-section for the heavy fermion production and LHC, as 
a function of the fermion mass. In section 5, we study the decay rates of these heavy 
fermions into Higgs and gauge bosons. We compare and contrast our model against the 
usual Type HI seesaw models with only one Higgs. We also show the consequences of 
non-trivial mixing of the heavy fermions on the fiavor structure of their decays. In section 
6, we discuss the decay rates and branching ratios of the Higgs decays. We probe issues 
on Higgs decays, which are specific and unique to our model. Section 7 is devoted to the 
discussion of displaced decay vertices as a result of the very long living in our model. In 
section 8, we list all possible final state particles and their corresponding collider signature 
channels which could be used to test our model. We calculate the effective cross-sections 
for all channels at LHC. We highlight some of the channels with very large effective cross- 
sections and discuss the standard model backgrounds. Finally, in section 9 we present 
our conchisions. Discussion of the scalar potential, the Higgs mass spectrum and the 
constraints from neutrino data on the Higgs sector is discussed in detail in Appendix A. 
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The lepton-Higgs couphng vertices are hsted in Appendix B.l, the lepton-gauge couphng 
vertices are hsted in Appendix B.2, and the quark-Higgs couphng vertices are hsted in 
Appendix B.3. 



2 Yukawa Couplings and Lepton Masses and Mixing 

We add three extra SU(2) triplet fermions to our standard model particle content. These 
fermions belong to the adjoint representation of SU(2) and are assigned hypercharge Y = 0. 
This makes each of them self conjugate. We will denote their Cartesian components a^ 




^: = , (3) 



c 

where i = 1, 2, 3 and \E'^ = \E''j . In the compact 2x2 notation they will be represented in 
our convention as 

s: = ^Esi-^, (4) 

where aj are the Pauli matrices. The right-handed component of this multiplet in the 2x2 
notation is then given by 



where 



^'i. = ^^^""^f^ E'l = E% (6) 

are the components of the triplet in the charge eigenbasis. The corresponding charge- 
conjugated multiplet will then be 



(5) 



^£ = Cmt = r^jr (7) 



The object which transforms as the left-handed component of the S multiplet can then be 
written as 



^Throughout this paper we denote particles in their weak eigenbasis by primed and mass eigenbasis by 
unprimed notation. 
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such that S',- = S' 



R^ 



^ Ri 



As discussed in the introduction, we include in our model a new SU(2) scalar doublet, 
$2, in addition to the usual standard model doublet $i. This new doublet couples only to 
the triplet fermions introduced above. The triplet fermions on the other hand are restricted 
to couple with only the new $2 doublet and not with $1. This can be ensured very easily 
by giving Z2 charge of —1 to the triplet fermions and the scalar doublet $2, and Z2 
charge +1 to all standard model particle^. The part of the Lagrangian responsible for the 
lepton masses can then be written as 



-A 



YiJ^Ml', + y^.M^R.L- + h.c. + ^M,, Tr 



J^'nt'l + h.c. 



(9) 



where L' and are the usual left-handed lepton doublet and right-handed charged leptons 



respectively, Yi and ¥■£ are the 3x3 Yukawa coupling matrices, and $2 
Higgs doublets $1 and $2 take Vacuum Expectation Value (VEV) 



($2) 



ia2^2- Once the 



(10) 



we generate the following neutrino mass matrix 




+ h.c, 



'111 



and the following charged lepton mass matrix 



Ri 



Ri 




(12) 
(13) 



Note that due to the imposed Z2 symmetry neutrino masses depend only on the new Higgs 
VEV v' while in the charged lepton mass matrix both the VEV's enter. The value of v' is 
determined by the scale of the neutrino masses and is independent of the mass scale of all 
other fermions. Therefore, the neutrino masses can be naturally light, without having to 
fine tune the Yukawas Is to unnaturally small values. 

The 6x6 neutrino matrix (ITT]) can be diagonalized to yield 3 light and 3 heavy Majorana 
neutrinos. The 6x6 unitary matrix U which accomplishes this is defined as 



T 







m 



T 
D 

M 



U 



Drr 






Dm 



and 




(14) 



•^We will break this Z2 symmetry mildly in the scalar potential. We discuss the phenomenological 
consequences of this Z2 symmetry and its breaking when we introduce the scalar potential and present 
the Higgs mass spectrum in Appendix A. 
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where mn = v 'Ys/V2, and 





( nil 













1712 







^ 


















Dm = 





Ms, 







I 








(15) 



Here rrii and Mj^- {i = 1, 2, 3) are the low and high energy mass eigenvalues of the Majorana 
neutrinos respectively. We reiterate that the primed and unprimed notations represent the 
weak and mass eigenbases respectively. The mixing matrix U can be parameterized as a 
product of two matrices 

U = W,U, (16) 

where Wiy is the matrix which brings the 6x6 neutrino matrix given by Eq. fill I) in its 
block diagonal form as 

^2-/0 m-^ \ / rh 



while Ui, diagonalizes fh^ and as 

rfrh 0\ (Dm 

U MJ^" ~ \ Dm 

The above parameterization therefore enables us to analytically estimate the mass eigen- 
values and the mixing matrix U in terms of W^y and Ui, by a two step process, by first 
calculating Wi, and then U^. Since the unitary matrix U has 6^ = 36 free parameters and 
the matrix 11^ has 2 x 3^ = 18 parameter, the matrix Wi, should have 36 — 18 = 18 free 
parameters. This matrix therefore can be parameterized as [29] 



where B = Bi + B2 + + ... and Bj ~ (l/ms)-', where is the scale of the heavy 
Majorana fermion mass. Using an expansion in l/m^ and keeping only terms second order 
or lower in 1/ms, we get 

" ~ [ -M-^TTiD 1 - iM-'mnm^M-'Y J ' ^ ' 

The light and heavy neutrino mass matrices obtained at this block diagonal stage are given 
by (only second order terms in l/mg are kept) 

fa^ = -ml,M-^mD, (21) 
M = M+- [mDm^^{M-^y + (M-^)*m^m|^) . (22) 
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Note that Eq. (!2T!) is the standard seesaw formula for the hght neutrino mass matrix, 
while Eq. ( !22|) gives the heavy neutrino mass matrix. These can be diagonalized by two 
3x3 unitary matrices Uo and t/s, respectively. This yields 

For the charged leptons we follow an identical method for determining the mass eigen- 
values and the mixing matrices. However, since the charged lepton mass matrix Mi given by 
Eq. f|T3l) is a Dirac mass matrix, one has to diagonalize it using a bi-unitary transformation 

^iSn. m)«-(? Z^J-^'- (^^) 

where m; = vYi, while Di and are diagonal matrices containing the light and heavy 
charged lepton mass eigenvalues. With the above definition for the diagonalization, the 
right-handed and left-handed weak and mass eigenbases for the charged leptons are related 
respectively as. 



Instead of using Eq. (124|) for the diagonalization, we will work with the matrices 

Ml Ml = SMIMi^ S\ and MiM} = T Mi^mI T\ (26) 

to obtain S and T respectively. As for the neutrinos, we parameterize 

S = WlUl, and T = WrUr, (27) 

where Wl and Wr are the unitary matrices which bring M^Mi and MiM} to their block 
diagonal forms, respectively, 

WlMlMiW, = ("^''p ~ \ , and W^mMWr = f ~ °, A . (28) 

^ ' \ Mh Mh I R t I R y Q MhMh j 

Using arguments similar to that used for the neutrino sector, and keeping terms up to 
second order in 1/ms, we obtain 

_ (l-mUM-^YM-^mn V2m^j,{M-^y \ , 

^ ^ -V2M-'mD 1 - M-^mDmUM~^y J ' ^ ' 

W - ( 1 ^mim'^^{M-yM-^\ 
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The square of the mass matrices for the hght and heavy charged leptons in the flavor basis 
obtained after block diagonalization by Wji and Wl are given by 



mim\ = mimj — 2mim^jy{M ^)*M ^mum], (31) 
+ mDm]mim^ij{M-^)*M-\ (32) 



and 



m\mi = mlmi — [rn}j;)M* ^M ^ m £)m\mi +h..c\ (33) 
M^H^H = M^M + M-^mDm^^M + M^mDm^^{M-^y + M-\mDm^^y{M~^y 

+ [M-\M-^yM-^mD{m]mi)m^^M-^M-^mDm^^M*-^M-^mDm^^M + h.cj;34) 

One can explicitly check that the masses of the heavy charged leptons obtained from Eqs. 
(!32|) and (IMI) are approximately the same as that obtained for the neutral heavy fermion 
using Eq. (122|) . Indeed a comparison of these equations show that at tree level, the 
difference between the neutral and charged heavy fermions are of the order of the neutrino 
mass and can be hence neglected. One-loop effects bring a small splitting between the 
masses of the heavy and neutral fermions of ~ 166 MeV. This allows the decay channel 
= S° + TT^ at colliders, as discussed in detail in [18, 19]. In this paper we will neglect 
this tiny difference and assume that the masses of all heavy fermions are the same. 

The matrices m\mi and M\jMh are diagonalized by Ui and t/^f giving, 

Ul={"' ".V (35) 



Similarly the rhirhj and MhM\j matrices are diagonalized by and and hence give, 

Ur={^: ur]- (36) 



VO Ui 

Finally, the low energy observed neutrino mass matrix is given by 

UpMNS = U}Uq. (37) 
Note that both Ui and Uq are unitary matrices and hence Upmns is unitary. 



3 A /i-T Symmetric Model 

As discussed in the introduction we wish to impose /i-r symmetry on our model in order to 
comply with the current neutrino data, which shows a preference for = and 623 = 7r/4. 
Henceforth, we impose the /x-r exchange symmetry on both the neutrino Yukawa matrix 
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Figure 1: Scatter plots showing variation of sin^ 612 (upper panels) and R = Am2i/|Am|]^| 
(lower panels) as a function of 04, an and a^. All Yukawa couplings apart from the 
one plotted on the x-axis, are allowed to vary freely. Only points which predict oscillation 
parameters within their current 3cr values are shown. Blue points are for mo = v'"^/ (2Mi) = 
0.95 eV while the green points are for rriQ = 0.006 eV. 
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Figure 2: Scatter plot showing the values of the Yukawa couplings which give all oscillation 
parameters within their current 3a allowed ranges. Allowed points are shown for mo = 0.96 
eV (blue), 0.006 eV (green) and 0.0021 eV (red). All Yukawa couplings apart from the 
ones plotted in the x-axis and y-axis are allowed to vary freely, in each panel. 



10 



0.001 



0.0001 



1e-05 



1e-06 



1e-07 



1e-08 




0.001 



0.0001 



1e-05 




1e-06 



Figure 3: Non-zero values of Ue3 and |0.5 — sin^^23| predicted when /i-r symmetry is 
broken. Shown are the oscillation parameters against the /i-r symmetry breaking parameter 
e = M3 — M2. Only points which reproduce the current neutrino observations within their 
3a C.L. are shown. The plot is generated at a fixed set of Yukawa couplings and heavy 
neutrino masses. 



Yy; and the Majorana mass matrix for the heavy fermions M. Therefore, the neutrino 
Yukawa matrix takes the form 



Yy 



^ CI4 (III (^11 

\a'ii as ag 



(3J 



In what follows we will assume (for simplicity) that a'n = an. Note that the /i-r symmetry 
does not impose this condition. It only imposes that Fs^j = Fsis Y-^^^ = Isg^. We 
have put a[^ = an in order to reduce the number of parameters in the theory. For the 
same reason, we assume all entries of Yj^ to be real. The heavy Majorana mass matrix is 
given by 



M 




(39) 



where without loosing generality we have chosen to work in a basis where M is real and 
diagonal. Here the condition M3 = M2 is imposed due to the ^ — r symmetry. 

The above choice of Yukawa and heavy fermion mass matrix lead to the following form 
of the light neutrino mass matrix 



,'2 



m 



Ml 



2 

n 

M2 



ail 



I a4 



(^u\ if: + 



M2 



"ll^^A/i Ah J 



Ml 



0.4 _|_ 0.6+ as 
Ml M2 

M2 

2a6ffl8 
M2 



ail 



Oj _|_ "6 +08 

Ml A/2 ) 



"11 I 
Ml ^ 



+ 



ril _|_ 

Ml ^ 

"11 I 
Ml 



2a6a8 
Af2 



Af2 



(40) 
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where we have used the seesaw formula given by Eq. (12T1) . which is vahd up to second 
order in l/m^- One can straightaway see from the above mass matrix that the scale of 
the neutrino masses emerges as ~ f'^a^/(2ms), where a is a typical value of the Yukawa 
coupling in Eq. (l38ll and the scale of heavy fermion masses. As discussed in the 
introduction, we restrict the heavy fermion masses to be less than 1 TeV in order that 
they can be produced at the LHC Therefore in principle, neutrino masses of ~ 0.1 eV 
could have been obtained with just the standard model Higgs doublet by reducing the 
Yukawa couplings to values ~ 10~^. However, this is usually considered as extreme fine 
tuning as there is no reason why the Yukawa couplings of the neutrinos should be so much 
suppressed, and the motivation for the seesaw mechanism is lost. In order to circumvent 
this, we introduced a different Higgs doublet $2, which couples only to the exotic fermions. 
On the other hand, Yukawa coupling of the standard Higgs $1 with the exotic fermions was 
forbidden in our model by the symmetry. Hence, only the VEV of this new Higgs doublet 
appears in Eq. (HO]) . Since this Higgs $2 is not coupled to any standard model particle, 
it could have a VEV which could be different. Therefore, we demand that f ' ~ 10^ eV in 
order to generate neutrino masses of ~ 0.1 eV keeping the Yukawa couplings ~ 1. We have 
checked that such low value of Higgs VEV is not in conflict with any experimental data. 
We will discuss in detail the scalar potential and the Higgs mass spectrum in Appendix A. 

We next turn to predictions of this model for the mass squared differences and the 
mixing angles. Since the neutrino mass matrix we obtained in Eq. ( HOl) has /i-r symmetry 
it follows that 

^13 = and ^23 = vr/4. (41) 

To find the mixing angle Q\2 and the mass squared differences Am^i and Amg^ 0, one needs 
to diagonalize the mass matrix m given in Eq. (HOj) . In fact, the form of rh in Eq. (HOl) is the 
standard form of the neutrino mass matrix with /i-r symmetry, and hence the expression 
of mixing angle 612 as well Am|^ and Am|]^ can be found in the literature (see for e.g. [30]). 
We show in Fig. [T]the variation of sin^6'i2 (upper panels) and R = Am2i/\Amli\ (lower 
panels) with the Yukawa couplings and ag. We do not show the corresponding 

dependence on since it looks almost identical to the panel corresponding to ag. The 
figure is produced assuming inverted mass hierarchy for the neutrino, i.e., Amg^ < 0. The 
neutrino mass matrix given by Eq. fHOl) could very easily yield Amg^ > and hence the 
normal mass hierarchy (see for e.g. [30]). However, for the sake of illustration, we will show 
results for only the inverted hierarchy in this paper. In every panel of Fig. [1], all Yukawa 
couplings apart from the one plotted on the x-axis, are allowed to vary freely. The points 
show the predicted values of sin^ 6*12 (upper panels) and R (lower panels) as a function 
of the Yukawa couplings for which all oscillation parameters are within their current 3a 
values [31], 

7.1 X IQ-^eV^ < Am^i < 8.3 x lO'^eV"^ , 2.0 x lO'^eV^ < lAm^J < 2.8 x lO'^eV"^ , (42) 

0.26 < sin^ 9^2 < 0.42 , sin^ 2^23 > 0.9 , sin^ < 0.05 . (43) 
''We define Am? — mf — m^. 
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For this figure we take Mi = M2 for simplicity and define mo = v'"^/ (2Mi). The blue points 
are for m = 0.95 eV while the green points are for tuq = 0.006 eV. 

Fig. [2] is a scatter plot showing the values of the Yukawa couplings which give all 
oscillation parameters within their current 3a allowed ranges given in Eqs. (H2ll and (H3il . 
Again as in the previous plot we assume Mi = M2, define mo = f'^/(2Mi) and show the 
allowed points for mo = 0.96 eV (blue), 0.006 eV (green) and 0.0021 eV (red). All Yukawa 
couplings apart from the ones shown in the x-axis and y-axis are allowed to vary freely, in 
each panel. 

Since the fi and r charged lepton masses are different, we phenomenologically choose to 
not impose the fx-r symmetry on the charged lepton mass matrij|§. Hence, without loosing 
generality, the charged lepton Yukawa matrix can be taken as 



Y, 



(Ye 

Y^ I , . (44) 

V 













/I 









1 




Ur ~ 





1 




Vo 









[0 








The masses of the light charged leptons can then be obtained from Eqs. fl3Tl) and/or fl33l) . 
For our choice of Yj^ and M, it turns out that me ~ vY^, , m^ ^ vY^, and 171^- ^ vYt-, if 
we neglect terms proportional to v''^. The mixing matrices Ui and Ur which diagonalize 
mj^m; (cf. Eq. fl55]) ) and rhirhj (cf. Eq. ( l3Ti) ) respectively, turn out to be unit matrices at 
leading order. 



(45) 



Finally, we show in Fig. [3] the impact of /i-r symmetry breaking on the low energy 
neutrino parameters. For the sake of illustration we choose a particular form for this 
breaking, by taking M3 7^ M2. Departure from /x-r symmetry results in departure of 
Ue3 from zero and sin^ 623 from 0.5. We show in Fig. [3] the f/es (left hand panel) and 
|0.5 — sin^ 0231 generated as a function of the symmetry breaking parameter e = M3 — M2. 
The figure is generated for 04 = —0.066, an = 0.171, Qq = 0.064, ag = 0.0037 and 
mo = 0.745 eV. We have fixed Mi = M2 = 299 GeV in this plot. For e = 0, /i-r symmetry 
is restored and both Ue3 and 0.5 — sin^6'23 go to zero. We show only points in this figure 
for which the current data can be explained within 3a. We note that for e > the curve 
extends to about M3 = M2 + 2.6 GeV, for this set of model parameters. For e < 0, the 
allowed range for e is far more restricted. 

We next turn our attention to the predictions of this model for the heavy fermion 
sector. The masses of the heavy fermions can be obtained using 1^ and M, as discussed 
in the previous section. The 6x6 mixing matrices, which govern the mixing of the heavy 
leptons with light ones, can also be obtained as discussed before. We will see in the next 
section that all the four 3x3 blocks of the matrices U, S and T are extremely important 



^We reiterate that our choice of the lepton masses and mixing are dictated solely by observations. 
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for phenomenology at the LHC We denote these 3x3 blocks as 



rr^fUn Uu\^f{W,)nUo {W,)i2U^\ .... 
U21 U22J \{W,)2iUo {W,)22Uj:J' ^""^ 



\T2X T22J \{WR)2lUr {Wn)22U^J' ^ 



The matrices W^, Wl and Wr have been given in Eqs. ( ]20l) . ( ]29i) and (!30ll respectively. 
These can be estimated for our choice of rriD, M and mi. In particular, we note that Su and 
Til are close to 1, while Uu is given almost by Upmns- The off-diagonal blocks f/12, U21, 
S12 and 5*21, are suppressed by ~ mjy/M, while T12 and T21 are suppressed by ~ m£)mi/M'^. 
Finally, the matrices ?722 = (^^1^)22^25 >S'22 = {Wl)22U^, and while T22 — f//f . To estimate 
these we need to evaluate first the matrices which diagonalize the heavy fermion mass 
matrices M, MjjMn, and MhM\j respectively. For M with /x-r symmetry, it turns out 
that 



t^." - f/.^ - I 72 -72 h (49) 
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thereby yielding 



U22 ^ \ ^ T2 ~T2\- (50) 
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S22 ^ --ij . (51) 



T22 ^ 71 • (52) 



This structure of f/^, Uf^ and f/;f (and hence of U22, S22 and T22), is an immediate conse- 
quence of the n-T symmetry in M and m^). This is an extremely new and crucial feature. 
To the best of our knowledge, this has not been pointed out in any of the previous Type 
III seesaw models studies. The main reason is that no study so-far has considered the 
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flavor aspect of Type III seesaw. As a result none of them considered imposing an un- 
derlying flavor symmetry group on the fermions such that the triplet fermion Majorana 
mass matrix and the Yukawa matrix would be /x-r symmetric. They assume that ?7s, Uj^ 
and are almost unit matrices since M is real and diagonal. However, this is true only 
if Ml <C M2 <S M3. Breaking of the yU-r symmetry either in or in M, will destroy 
this non-trivial form for [/s? and Uj^. But having /i-r symmetry in both and M is 
both theoretically as well as phenomenologically well motivated. We will see later that this 
non-trivial form of the matrices U^,-, Uj^ and f/^f will lead to certain distinctive signatures 
at LHC. 

One should note that while f/^, Uj^ and Uj^ have the form given by Eq. ( H9l) . Ui and 
Ur ~ J, though both M and Yi were taken as real and diagonal. The main reason for this 
drastic difference is the following. While we take exact /i-r symmetry for M, for Yj we 
take a large difference between and Y^ values. This choice was dictated by the observed 
charged lepton masses. 

Finally, a comment regarding the extent of deviation from unitarity in our model is 
in order. It is clear from the discussion of the previous section and Eq. fl20|) that the 
deviation from unitarity of the light neutrino mixing matrix is oc m|,/m| ~ rrij^/ms, where 
m,^ and ms are the light neutrino and heavy lepton mass scales respectively. Therefore, 
an important difference between our model and the usual GUT Type III seesaw models 
is that the extent of non-unitarity for our model is much larger. This will result in larger 
lepton flavor violating decays in our case. Detailed calculations for lepton flavor violating 
radiative as well as tree level decays of a generic Type III seesaw model have been published 
in [14, 15] and we do not repeat them here. The authors of these papers have also worked 
out the current constraints on the deviation from unitarity. One can check that even for 
100-1000 GeV mass range heavy leptons, the predicted non-unitarity and lepton flavor 
violating decay rates in our model are far below the current experimental bounds. In fact, 
the predicted decay rates can be seen to be far below the sensitivity reach of all forthcoming 
experiments. 

4 Heavy Fermion Production at LHC 

The triplet fermions couple to the standard model particles through the Yukawa couplings 
as well as gauge couplings. We give in Appendix B, the detailed Yukawa and gauge 
couplings of the neutral and charged heavy fermions with the standard model leptons, 
vector bosons, and Higgs particles. We have kept the masses of the heavy fermions in the 
100 GeV to 1 TeV range. Therefore, it should be possible to produce these fermions at 
LHC. In this section, we will study in depth the production and detection possibihties of 
the heavy leptons in our Type HI seesaw model. Compared to the earlier papers, there 
are two distinctly new aspects in our analysis - (i) presence of two Higgs doublet instead 
of one, leading to a far more rich collider phenomenology, (ii) presence of /i-r symmetry in 
our model. 

The heavy triplet fermion production at LHC has been discussed by many earlier papers 
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[17-21]. At LHC we will be looking at the following production channels 

The exotic fermions have both Yukawa couplings to Higgs as well as gauge couplings 
to vector bosons. Therefore, they could be in principle produced through either gauge 
mediated partonic processes (left diagram) or through Higgs mediated partonic processes 
(right diagram) 

q/q' S± q/q' E± 




q S7 g TP/ 

However, it turns out that the vertex factors for the couplings of heavy fermions to gauge 
bosons which are relevant for the formers production, viz., T'^YT Zj^ and Y^Ys^W'^ , are 
much larger than those involving the Higgs, viz., Y'^Y~ /h^ / and Y^Y'^H'^ . To illus- 
trate this with a specific example, we compare the Y'^Y^ Z coupling given in Eqs. (]B16|) 
and (IB 171) with the Y'^Y^h^ coupling given in Tabled It is easy to see that the 
coupling has terms proportional to T22T22 and S'225'22, while the Y'^Y~h^ coupling depends 
on terms which have an off-diagonal mixing matrix block. Since the off-diagonal mixing 
matrix blocks are much smaller than the diagonal ones (as discussed in section 3), it is not 
surprising that the couplings of two exotic fermions to the Higgs particles are much smaller 
than to the gauge bosons. Hence the heavy exotic fermions will be produced predominantly 
via the gauge boson mediated processes. For the heavy fermion production cross-section at 
the LHC, we chose CTEQ6L [32] as the parton distribution function (PDF) and partonic 
center of mass energy as the renormalization and factorization scale. We have explicitly 
checked that the production cross-sections do not change much with the PDF and scale. 
All cross-sections in this paper are calculated using the Calchep package [33]. 

In Fig. Hlwe show the production cross-section for Y~Y^ (bold dotted line), (solid 
line), and S+E" (fine dotted line), at LHC as a function of the corresponding heavy fermion 
mass. It is straightforward to see that the (and E'°S'V±) couphngs are absent. A 

very small Y^Y^Z is generated through mixing from the u'^u'^Z coupling. However, this is 
extremely small. Hence, production through gauge interactions is heavily suppressed 
and is not shown in Fig. |H One can see that the production cross-sections of the heavy 
fermions fall sharply with their mass. More precisely, (t(S^S^) = 112 fb, a{Y'^Y'^) = 206 
fb and a{Y-Y^) = 95 fb, for Ms, ~ 300 GeV. However, for Ms, ~ 600 GeV this quickly 
falls to a(E±E^) = 6 fb, a{Y+Y'^) = 13 fb, and a{Y-Y°) = 4 fb. Therefore, it is obvious 
that the lightest amongst the three generation of heavy fermions will be predominantly 
produced at the collider, and will hence dominate the phenomenology. One can check that 
the production cross-sections that we get is almost identical to that obtained in earlier 
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Figure 4: Variation of production cross section of S^*", S° with the mass of exotic leptons. 



papers [19,20]. This is not unexpected since our model is different from all the earlier 
models in the Higgs sector. However as discussed above, it is the gauge interactions which 
predominantly produce the exotic leptons. The gauge-lepton couplings in our model is 
same as in the earlier works. Since the heavy fermion production comes mostly from the 
gauge mediated sub-processes, we get the same production cross-sections as other papers. 

5 Heavy Fermion Decays 

Once produced at LHC, the heavy fermions will decay to different lighter states due to its 
interaction with different standard model particles. In particular, they could decay into 
light leptons and Higgs due to their Yukawa couplings, or to hght leptons and vector bosons 
due to their gauge interactions. The light leptons could be either the charged leptons or 
the neutrinos. The Higgs could be either the neutral Higgs h^, H^, A^, or the charged 
Higgs H"^. The gauge bosons could be either or Z. The Higgs and gauge bosons 
would eventually give the final state particles in the detector, which will be tagged at the 
experiment. This will be studied in detail in the following sections. Here we concentrate 
on only the two body decay rates and branching ratios of the exotic leptons and S''. 
AU possible vertices and the corresponding vertex factors for the Yukawa interactions of 
and are given in Tables [9|, [TOl [TTl The vertices and vertex factors for the charged 
and neutral current gauge interactions can be found in Appendix B.2. Presence of two 
Higgs doublets and fi-r symmetry in 1^ and M will have direct implications in the partial 
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decay widths for different decay processes. 



5.1 Decay to Light Leptons and Higgs 

In this subsection, we perform a detailed study of all two-body decays of these fermions 
into a lepton and a Higgs. Since we have two Higgs doublets in our model, we have a pair 
of charged Higgs if^, and three neutral Higgs - and are CP even, while is CP 
odd. The Higgs mass spectrum and mixing is given in Appendix A. The decay width F for 
Sj — > IjX is given by 



32n 



My 



X A 



■311 



(53) 



where M^^, Mx and rrti. are the masses of S°/Sf , X and /j, respectively, where X is the 
relevant Higgs involved. The Ij could be either a charged lepton or a neutrino. For the 
charged Higgs B.^ mode, and the neutral Higgs hb^ and mode, the factor Aji is given as 



A 



IS 



2 I |/^X,_R< 



(m; 



Mi 



) M, 



(54) 



while for the CP-odd neutral Higgs A^ the factor is 



= (i(Q^'%r+i(QE'%r)(i- 



(Ml - Mf 



Ml 



■iX,L\ 



Mv: 



(55) 



In the above equations {C^^)/{C^") are the relevant vertex factors given in Table [9| [TO] 
anddH and i, j represents the generation. In all numerical results that follow we will fix the 
model parameters (Yukawa couplings and entries of M mass matrix) to their values given in 
Tabled This set of model parameters yield Am|i = 7.67x 10"^ eV^ Am|i = -2.435 x 10~^ 
eV^ and sin^ 6'i2 = 0.33. Of course 6^13 = and 623 = ir/A. Throughout the rest of the 
paper, we also take the value of M^o = 40 GeV, Mho = 150 GeV, M^o = 140 GeV and 
Mh± = 170 GeV. Also, for all cases where we present results for fixed values of the heavy 
fermion masses, we take M^^ = 300 GeV and Ms2 = M^g = 600 GeV. 



5.1.1 ^l^hyH'^/A^ 

Let us begin with the decay of heavy charged leptons into light charged leptons and neutral 
Higgs. The Higgs concerned in this case could be h^, H^, or A'^. We start by probing the 
decay rate Sf — > ifh^- From Eq. ([S3D we see that the decay rate is governed by the factor 
Aji, which in turn depends on the vertex factors given in Tabled The vertex factors are 
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04 


ae 


as 


an 


mo/eV 


^2 


Ms 
Ml 


0.145 


0.097 


0.109 


4 X 10-^ 


2.356 


2.0 


2.0 



Table 1: Model parameters used for all numerical results in section 5 and 6. This set 
of model parameters yield Am2i = 7.67 x 10"^ eV^, Am|^ = —2.435 x 10~^ eV^ and 
sin2 6'i2 = 0.33. Of course 6113 = and 6^23 = 7r/4. Parameter mo = t;'V(2Mi). 
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Figure 5: Variation of r(Sj — > / • h^) with M^. 



given in terms of the 3x3 block matrices Sab and Tab, where a,b = 1, 2. We have seen in the 
earher sections that S12, T12 and T21 are heavily suppressed - the first one by 0{mD/M) 
and T12 and T21 by 0{{mimD)/M^). The vertex factors also depend on the Higgs mixing 
angle a. In Appendix A, we have shown how the neutrino mass constrains the neutral 
Higgs mixing such that sin a ~ 10~^ and cos a ~ 1. Therefore, for the Sf — >■ Z^/i*^ decay 
the dominating vertex factor is 

Cf4i ^ -^SIy^T22 cos a. (56) 
v2 

We have seen in Eq. f H7|) that ~ 1 if we neglect terms of the order of 0{m'jy/M^). 
Therefore, 



/ 04 V^cbii 

an 75(^6 + as) ^(as-cte)!. (57) 
\aii ^(ae + as) -^{ae - as) 
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Figure 6: Variation of r(Sj —>■ IjH) with M-£. 



From Eq. fl57p we can see that {Ci± E±)i3 — 0. In fact one can check that this happens 
because T22 given by Eq. (1521) has a specific form, which is due to fi-r symmetry. The 
consequence of this is that decay of S3 — e~h^ will be forbidden to leading order. Also 
note from Eq. fl57j) that the decay rate of all heavy charged fermions into fi^ is predicted 
to be exactly equal to their decay rate into r^. This is also an obvious consequence of the 
fi-T symmetry. 

The partial decay widths for S~ l~h^ from an exact numerical calculation in shown 
in Fig. [5], as a function of the heavy charged fermion mass. The thin blue lines are decay 
into e~, while the thick green lines are for decay into /i~/r~. As expected, we notice the 
following two consequences of /i-r symmetry 

• Decay rate of S3 e^h^ is almost zero. 

• The decay rate of the heavy fermions into fi~ is exactly equal to that into r^. 

We can also see that for Sj" decay, the decay rate into e~ is about 5 orders of magnitude 
larger than into fi^/r^. The trend is opposite for Sj decay, where the decay is predomi- 
nantly into fi^ /t^ . Both of these features can be understood from Eq. fl57|) and the values 
of the Yukawa couplings taken (cf. Tabled]). S^f decay into e~ and fi~ /t~ is proportional 
to al and af^, respectively. The ratio of the decay widths seen in the figure matches the 
ratio 04/(^11) ~ 10^. Similarly, one can check that for Sg decay, the corresponding ratio 
is 4:a1i/ (ag + ag)^, which agrees with the middle panel of Fig. O Finally, the fact that the 
decay rate of S3 fi~h^ is less than that of S2 fi~h^ can also be understood in terms 
of Eq. (l57|l and the Yukawa coupling values taken for the calculation. 
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Figure 7: Variation of r(S° ^ ljH+) with Ms, 



We next turn to the decay width for Ef — > ifH^. Expression for the decay rate is same 
as that given by Eq. (!53l) except that now M/jO is replaced by the mass M^o. For this 
decay channel the Aji factor is dominantly given by 

C&^^^SlYiT,2sma. (58) 

Note that compared to the effective vertex factor for l^h^ given in Eq. (jMD, 

the effective vertex factor given above for Sf — > ifH^ is suppressed by sin a. Since 
sin a ~ 10^^, the decay rate of Sf into are heavily suppressed. We show in Fig. [6] this 
decay rate calculated from exact numerical results. Comparing Fig. [5] with Fig. [6], we see 
that decays into are suppressed by a factor of about ~ 10^°, as expected from the order 
of magnitude estimate. Therefore, we can neglect Ef ^f^^ practical purposes. 

From Table [H] it is easy to see that the decay rate Sf ^f^^ be almost identical 
to that that predicted for Sf — > ifh^- The vertex factors for the two process are the same 
and hence the only difference could come from the difference between the Higgs masses. 
However, it is easy to see from Eq. fl53l) that the effect of the Higgs mass on the decay rate 
is not very significant, especially for relatively heavy fermions. 
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5.1.2 TP^l^H^ 



The decay rate for this channel is also given by Eq. (!53l) . and is governed primarily by the 
vertex factor 

Cj^i^ ^ ^SlY^U;, cos p. (59) 

Recall that cos/3 ~ 1. As discussed before, the matrix U22 displays features similar to the 
matrix T22. Therefore, the form of dominant vertex factor for this case is similar to that 
given in Eq. (!57j) . The corresponding decay rates are shown in Fig. [7l All features seen for 
Sj" —>■ Ijh^ is also seen here. Decay channel Eg — > ejH^ is forbidden. Decay rates to fi^ is 
equal to decay rate to t"^. The huge hierarchy in the decay rates of E5 and E2 into e and 
fi/r are also present due to same reason as given for S~ — > I'h'^ decays. The decay rate 
and flavor structure for the final state charged leptons is therefore seen to be same here 
as for the decay of charged heavy fermions into charged light leptons and h^. However, in 
this case we have a charged Higgs in the final state and it should be easy to tag this and 
differentiate the two processes in the detector at LHC We will also discuss in the following 
sections that the also has a much longer lifetime than H^, which can be observed in 
the detector. In addition, the heavy lepton itself is charged in one case and uncharged in 
the other. The two processes should hence be separable at the collider experiment. 



5.1.3 z/ /lO/iyo/AO 

We next turn to the decay channels with a light neutrino in the final state. This will give 
missing energy in the final state. Decay of the neutral S° will create a neutrino and a 
neutral Higgs. As in the case of decay of to charged leptons and neutral Higgs, one can 
check from Table [10] that the decay to the Higgs is heavily suppressed due to smallness 
of U21 as well as the sin a term. However, decay to /i° is driven by the vertex factor 

cf,^o'' = lulY^U;,cosa. (60) 

For the decay S° ^3-^ we find from Table [TOl that the dominant vertex factor is 

Cti^ = --^ly^Vl2^o^i^- (61) 

Since cos/5 ~ cos a, the decay rate and flavor structure for this channel will be similar to 
what we found for the S° z/j/i" channel. The main difference comes in the difference 
between the masses of the /i" and Higgs. 

For the ^jH^ decay, one can see from Table [TO] that the vertex factors for both as 

well as Pr vertices, are suppressed by sin a. Therefore, this decay rate can be neglected. 
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5.1.4 E^^uH^ 



From Table [TT] the vertex factor for this decay will be 

Cf^l'^ ^ f/^ilf ^22 COS (3. (62) 

As we have seen in section 3, the structure of 5*22 is very similar to that of 1/22- Hence, 
a comparison of the vertex factor for this process with the one from E° Ujh^ and 
S° — i> i^jA^ shows that all three will have decay rates of comparable magnitude, modulo 
the difference in the masses of the scalars h^, and H^. Since we assume masses of 

and as 40 GeV, 140 GeV and 170 GeV respective, the decay rate for uH^ 
is predicted to be the lowest. 

5.2 Decay to Light Leptons and Vector Bosons 

The exotic heavy leptons have gauge interactions. Therefore, it is expected that they will 
also decay into final state particles with vector bosons, and Z. The decay width F for 
Sf if / and — > uV in the mi = limit is given by 

(63) 

where Cj±^ and C^±^ are the relevant vertex factors given in Appendix B.2, and My is the 
mass of the vector boson involved. The dominant vertex factor relevant for l^Z and 

S*^ — > l^W^ in terms of Y^,, M, v' and the mixing matrices are given respectively by 

Cf4^ c^^^ uMM-'Ut, and C^^l^ ^ g Ui^Y^M-'U^. (64) 

For the other two channels TP —>■ vZ and Ys^ —>■ vW^ , they are given respectively by 

C',i^:^{9c^ + 9's^)uUiNsY^M-'U,,, and C,^!''' ^ ^7 f/jMiV5^E^M-if/,^.(65) 

As mentioned before, the gauge interaction part of our model is identical to that for the 
one Higgs doublet Type III seesaw considered earlier. Some of these vertex factor^ can 
therefore can be seen to agree with that given in [15]. The only difference is that we 
include the matrices Ui, and f/s in our general expressions, while these were taken as 
unit matrices in [15]. 

5.3 Comparing E^/^ Decays to Higgs and Gauge Bosons 

In Fig. Owe show the decay rates T,i uiW^ (long-dashed blue line), Ti Z (dot- 

dashed green line), ^ viZ (dot-dashed magenta line), S'j' e~W^ (thin solid red line), 

^Vertex factor for YP vZ is not given in [15]. 
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Figure 8: Comparison of the decay rate of the heavy fermion into (i) Higgs and (ii) vector 
bosons, in our two Higgs doublet model. 



—>■ e~H^ (dotted maroon line), Sj" viH~ (dashed violet line), and e~h^ (thick 

solid dark green line). One can clearly see that all decays to gauge bosons are suppressed 
with respect to decays to (and A^) and by a factor of more than 10^°. The reason 
for this can be seen by comparing the vertex factors involved in decays to Higgs and 
(cf. Eqs. (156]), ([S9l), (feni), ([61]), ([62])), with decays to gauge bosons (cf. Eqs. ^ and 
(I65p ). It is clear that while the former vertex factors do not have any suppression factor, 
the latter are all suppressed by f'/M^. Another important difference between the decay 
rates to Higgs given in Eq. fl53l) . and gauge bosons given in Eq. f l63]) . is in the kinematic 
factors. Comparison of the two equations reveals that (for = 0), there is an additional 
factor of (2 + M|/My.) for the gauge boson decays. This extra M| in the numerator cancels 
out the 1/M| in the denominator coming from the square of the vertex factors. However, 
the suppression of the gauge boson decay rates due to g^v''^ jMy oc v''^ jV'^ remains, where 
V^2 _ _|_ gince we have taken v' ~ 10~^-10~'' GeV, the decays to gauge bosons 
are suppressed by a factor of ~ 10^°-10^^ compared to the decays to Higgs. Therefore, 
branching ratios of the heavy fermion decay to and Z can be neglected in our model 
and we concentrate on only decays to Al* and in our next section. Note that the 
decay to is also suppressed by a factor of 10^°-10^^, as was also pointed out earher. 
We had seen that this suppression is due to sin^ a coming from the vertex factor for this 
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process. Since sin^ a ~ 10 we find that the decay rate for this case is of the same order 
of magnitude as the decays to the gauge bosons. Hence, this is also neglected henceforth. 

5.4 Comparison Between One and Two Higgs Doublet Models 

It is pertinent to compare the two-body decays of the heavy fermions in our two Higgs 
doublet model with the usual Type HI seesaw models considered earlier which have one 
Higgs doublet. The expressions for heavy fermion decays to Higgs and gauge bosons in 
the one Higgs doublet models have been given before in the literature [18-21], and we give 
them here for the sake of comparison. The decay rates to gauge bosons in the one Higgs 
doublet model is given as (for rrii = 0) 

,HDM^j.o^^Z) ^ ^!i^(l-il|)2(i + 2i^|), (66) 

^IHDM^^O _^lT^r±^ ^ ^!i^(l_i^)2(l + 2^), (67) 

^ ^ 327r ^ My^' ^ Me^^' ^ ' 



^,nuM^^±^l±2) ^ ^!i^(l-il|)2(i + 2i%), (6J 



^.mmf ^^.i ^ ^^±^ ^ :^^(l_i!l^)2(i + 2^). (69) 
^ ^ 167r ^ Ms^^ ^ My^^ ^ ^ 

where A is the triplet fermion - lepton doublet - Higgs doublet Yukawa coupling in the one 
Higgs doublet model, and all mixing terms are neglected. This should be compared with 
the corresponding expression given in Eq. (1631) . which on neglecting all mixing and hence 
flavor effects reduces to (for mi = 0) 

p2™/(^0^^^) ^ ^li^!^(l_il|)2(i + 2^), (70) 

^2HDM^^, ^l^^±^ ^ 4^^(1-^)^(1 + 2^), (71) 

Ms Ms 
Ms 



1/1 Ms 


/2 

V 


647r 




1/1 Ms 


V 


327r 




1/1 Ms 


/2 

V 


327r 




1/1 Ms 


l2 

V 


167r 





,4- Vv / IVly .0, 



(73) 



where V"^ = v"^ + v'"^ is the electroweak breaking scale. The scale of the Yukawa coupling 
constants and VEVs are fixed by the neutrino mass my ~ A^K^/Ms for the one Higgs 
doublet model and m^ ~ ?/|t>'^/Ms. If one replaces A^ and y\v''^ jV"^ with myM^jV'^ in 
both set of expressions, one can see that the the decay rates of heavy fermions into gauge 
bosons are identical for both models. 

The rates for decay into Higgs for the one Higgs doublet model neglecting flavor effects, is 
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given by (for rrii = 0) 



.IHDAU^O ,^^0) ^ (74) 

^ 647r ^ My^^ ^ ^ 



ri™^(E± ^ /^i/O) ^ (75) 

327r Ms 

For the two Higgs doublet model, the corresponding decay rates are given by Eq. (1551) . 
which on neglecting all flavor effects reduces to (for mi = 0) 

r2HZ.M(^o _ ^^0/^0) ^ y| cos^aM, jg^ 
r2HDM(^± ^ ^±^0/^0) ^ cos^-aiw. ^^_^^,^ ^^^^ 





cos^ a Ms 




647r 


yl 


cos^ a Ms 




327r 


1/1 


sin^ a Ms 




647r 




sin^ a Ms 


327r 



M 



s 



where the first two expressions are for decays to or and the last two for decays to H^. 
Again, for the same value of Ms ~ 100 GeV in both models, one requires A ~ 10~^-10~^ 
for the one Higgs doublet model in order to produce rrii, ~ 0.1 eV, while ?/s ~ 1 for our 
two Higgs doublet model. Therefore, clearly 



^2HDM^^± ^ ^±/,0/^0) _ ^Qll ^ pmm/(^± ^ ;±^0)_ 

Hence, the the exotic fermions decay about 10^^ times faster in our model compared to 
the one Higgs doublet models]. This could lead to observational consequences at LHC. In 
particular, authors of [19] talk about using "displaced vertices" as a signature of the Type 
III seesaw mechanism. In our model the lifetime of the exotic fermions is a factor of 10^^ 
shorter and so will be the gap between their primary production vertex and the decay 
vertex. Our model therefore predicts no displaced vertex for the heavy fermion decays. In 
addition, decay to are predominant. The decay predominantly into hh pairs, but with 
a very long lifetime, as we will discuss in section 6. This will give a distinctive signature 
of our model at LHC. We will discuss displaced vertices from decay in section 7. 



5.5 Flavor Structure and the Decay Branching Ratios 

Finally, we present the branching fractions of the heavy fermion decays. Table [2] shows the 
branching fractions for the E^, while Table [3] gives the branching fraction for TP decays. 

^We reiterate that the decays YP — > vH'^ and YP — > l^H^ are suppressed by the sin^ a ^ 10^^^ factor 
and hence turn out to be comparable to the decay rates in the one Higgs doublet model. However, the 
branching ratio to this mode is negligible and can be neglected. 
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Decay modes 




^2 


^3 




0.363 


0.473 


0.473 




0.247 


2.28 X 10-6 


0.0 


/i± AO 


2.3 X 10-6 


0.125 


0.125 


r± AO 


2.3 X 10-6 


0.125 


0.125 




0.389 


2.5 X 10-6 


0.0 




3.6 X 10-6 


0.139 


0.139 


r± /lO 


3.6 X 10-6 


0.139 


0.139 



Table 2: Decay branching fractions of S^, and Sf for Mft0=40, Mj:/o=150, Mh± = 170 
GeV and M40 = 140 GeV. We have taken model parameters Mi = 300 GeV and M2 = 
M3 = 600 GeV. 



Decay modes 


SO 


so 


so 
^3 




0.368 


4.3 X 10-6 


0.0 




3.4 X 10-^ 


0.236 


0.236 




3.4 X 10-s 


0.236 


0.236 




0.243 


0.250 


0.250 


Z//lO 


0.386 


0.277 


0.277 



Table 3: Decay branching fractions of S?, SO and SO for M;,o=40, Mho=150, M^± = 170 
GeV and = 140 GeV. We have taken model parameters Mi = 300 GeV and M2 = 
M3 = 600 GeV. 



For the channels with neutrino in the final state, we give the sum of the branching fraction 
into all the three generations, as observationally it will be impossible to see the neutrino 
generations at LHC. We do not show decays to gauge bosons and as they are suppressed 
by a factor of 10^^ with respect to the decays into W ^ Al^ and . As a result of the inherent 
/i-r symmetry in the model, S3 '''' decays to electrons is strictly forbidden and branching 
ratios of their decay into \i and r leptons are equal. We find that due to the form of f/22, 
5*22 and T22 given in Eqs. (!50|) . (IST!) and (!52|) . S^*^" decays to electrons is also negligible and 
their probability to decay into /i and r leptons is equal. We also find that the branching 
fractions of Sf is equal to the branching fractions of Sf , and similarly for the neutral 
heavy fermions. We also notice that S^''° decays only to electrons and their decay to \i 
and r lepton is almost zero. This as pointed out before, comes due to the constraint on the 
Yukawa couplings from the low energy neutrino oscillation data. The difference between 
the branching fraction to /^o, Al^ and is mainly driven by the difference in the masses 
which we have chosen for these Higgses. 
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6 Higgs Decay 



In the previous section we concluded that the heavy fermions will all decay into 
or H"^. We next turn to the subsequent decay of these Higgs particles. We concentrate 
on the possible decay modes of h'^, and and tabulate only those few which have 
significant branching ratios. The branching ratios obviously depend on our choice for the 
Higgs masses as well as our choice of the mixing angles a and (3, which appear in the 
coupling. The part of the Lagrangian containing the interaction terms of the Higgs with 
the leptons and quarks are given in Appendix B. The interaction of Higgs fields with 
the gauge fields comes from the Higgs kinetic terms and is the same as the general two 
Higgs doublet model. Possible decay channels for the charged Higgs involve the and 
the neutral CP even Higgs. It is well known that in the two Higgs doublet model, the 
— — coupling is proportional to sin(/5 — a), whereas — — coupling 
is proportional to cos(/5 — a) [23]. In Appendix A, we have shown how constraint from 
neutrino mass drives sin a ~ sin (3 ~ 10~^. Therefore, in our model — > W^H^ is always 
suppressed, irrespective of the Higgs mas^. In fact, the only decay channel possible for 
the charged Higgs in our model is H"^ W^h^, for which the decay branching fraction 



BR{H^ W^h^) = 1.0. (80) 

The next decay into either qq' pairs or l^ui/ui pairs with the following decay branching 
fractions 

BR{W^ qq') = 0.67, 
BR{W^ ^ e^Ue/ue) = 0.11, 
BR{W^ ^ fx^uju^,) = 0.11, 

BR{W^ T^Vr/Vr) = 0.11. (81) 

The branching fractions of the neutral Higgs and A^ are given in Table HI Though 

is almost never produced through heavy fermion decays in our model, we have included 
them in the table for completeness. We find that the neutral Higgs decay to bb pairs almost 
90% of the times. The second largest decay fraction is to rr pairs, while decays to cc 
happens less than few percent of the times. In our following sections where we look for 
collider signatures, we will consider (and A^) decays to only bb and rf pairs. 

Finally, a short discussion on direct production of h^, without involving the heavy 
fermion decays, is in order. In our model, the lightest Higgs has a mass as low as 40 GeV. 
This might appear to be a cause of concern, given that such a Higgs was not observed at 
LEP. However, it is easy to see that this Higgs mass is not excluded by the direct Higgs 
searches at LEP-2. This is because the coupling corresponding to Z — Z — vertex is 
given by {gMz/ cos 9^) sin(/5 — a). Since in our model sin(/5 — a) is almost zero, the LEP-2 
bound on Higgs mass does not pose any serious threat to our model, irrespective of the 
mass of h^. 



^One can see that this channel is also kincmatically forbidden for our choice of the Higgs masses, 
whereby with Mh± — 170, it is impossible to create an on-shell pair oiW^ and H°. 
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Decay modes 






A" 


bb 


0.89 


0.87 


0.87 


Tf 


0.07 


0.09 


0.09 


CC 


0.04 


0.04 


0.04 



Table 4: Decay branching fractions of and for = 40 GeV, Mho = 150 GeV, 

and Mao = 140 GeV. 

7 Displaced Decay Vertex 

Amongst the most significant difference of our model with the usual Type III seesaw 
model are the decay lifetimes of the heavy fermions and hP (as well as A^). We had seen in 
section 5 that the total decay rate for 300 GeV is about 4 x 10~^ GeV. This gives the 
corresponding rest frame lifetime as 4.97 x 10~^^ cm. The lifetime for is similar. One 
can check that for the usual one Higgs doublet models, the rest frame lifetime for the heavy 
fermions is ~ 0.5 cm for rUi, — 0.1 eV and Ms ~ 100 GeV, which is rather large. The 
authors of [19] therefore proposed that the displaced decay vertex of heavy fermion could 
be a typical signature of the one Higgs Type III seesaw model. Clearly, for our model 
with two Higgs doublets, the decay lifetime is 10^^ times smaller and hence we predict 
no displaced vertex for the heavy fermion decay. This can be used as a distinguishing 
signature between the two models. 

Another very important and unique feature of our model is the very long lifetime of our 
neutral Higgs /i", which comes due to the smallness of sin a. In fact, since sin a ~ 10~^, the 
lifetime for in our model is 10^^ times larger compared to the standard model Higgs. In 
particular, the total decay rate is 4 x 10~^^ GeV. This gives a rest frame lifetime of 
4.97 cm. For a hP with 200 GeV of energy, the lifetime in the lab frame is seen to be close 
to 25 cm. Therefore, we expect a big gap between the decay vertices of the heavy fermion 
and the This displaced /i° decay vertex should be detectable at the LHC detectors 
ATLAS and CMS. 

We would like to make just a few qualitative remarks about the prospects of detecting 
the displaced /i° vertex. Like stressed many times before, the /i° decay predominantly into 
bb pairs. While ^-tagging is a very important and standard tool for coUider experiments, 
and while both ATLAS [34] and CMS [35] have been developing algorithms for tagging 
the b, there is an additional complication with 6-tagging in our model which should be 
pointed out here. Since the /i° lifetime is a few 10s of cm in the lab frame, it is expected 
to decay inside the sihcon tracker of ATLAS and CMS. In particular, the pixel tracker of 
CMS and ATLAS which are only few cm from the center of the beam pipe, will miss the 

decay vertex. However, the silicon strip trackers would be useful in observing the 6-jets. 
The tracks from the primary and secondary vertices of the 6-hadron should be seen. In 
addition, one could use the two other standard tools for tagging the 6-jets. Firstly, one 
could the tag the lepton in the jet coming from the semi-leptonic decays of the 6-hadron. 
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These leptons are expected to have smaller px compared to the ones coming from W and 
Z decays, and hence this is called soft-lepton tagging [34,35]. More importantly, one could 
construct the invariant mass distribution of the 2 6-jets. This should give us a sharp peak 
corresponding to the h^. We therefore expect that ATLAS and CMS should be able to 
detect the displaced hP decay vertex. This would give a characteristic and unambiguous 
signal of our model. 

Note that while the lifetime of is constrained to be large due to the smallness of sin a 
alone in our model, things are slightly more complicated for the lifetime of A^. This is 
because in principle could decay through the mode A^ — > Zh^. While this is forbidden 
kinematically for the ^4° mass we assume, one could argue that for a large enough mass 
for the hfetime of ^4*^ could be shorter. However, we stress that even if A^ decays fast 
into hP, that would still produce a displaced vertex, since the hP would still have a very 
long lifetime. 

8 Model Signatures at the LHC 

Having discussed in details the production and subsequent decays of the exotic fermions, 
as well as the decay branching fractions of the intermediate Higgs into final state particles, 
we next describe the signatures of the two Higgs doublet Type HI seesaw model at the 
LHC. We will present a comprehensive list of final state particles and their corresponding 
collider signatures. 

The most important characteristics of our model are the following: 

1. Presence of /i-r symmetry in and M. This is expected to show-up in the flavor 
of the flnal state lepton coming directly from the E^/'^ decay vertex. 

2. Presence of two CP even neutral Higgses and H^), one CP odd neutral Higgs 
(^4°), and a pair of charged Higgs (H^). 

3. Predominant decay of the heavy fermions into hght leptons, and A° or H"^. 
Decays into and gauge bosons almost never happen. 

4. Very short lifetime for the heavy fermion due to the very large Yukawa couplings. 

5. Predominant decay of and A^ into bb pairs 89% and 87% of the time, respectively. 
They decay also into rf 7-9% of the time. 

6. Very large lifetime of /i° and A^. 

7. The Higgs decays into W^hP and almost never into W^H°. 

8. Short predicted lifetime for H^. 

In what follows, we will use these model characteristics to identify distinctive final state 
channels at the collider. We identify all possible channels in the collider for our model 
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and calculate the respective effective cross-sections. The results are given in Tables [5l [6] 
and [71 We will also discuss some of the most important channels and the characteristic 
backgrounds, if any, associated with them. In this section we have only given results for 
effective cross- sect ions for the decay of E^''° with M^^ = 300 GeV. Results for the other 
heavy fermion generations can be similarly obtained. 

8.1 Signatures from S+S decays 

We give in Table [5] all possible collider signatures coming from the decay of S+S^ pairs, 
for our two Higgs doublet Type III seesaw model. In the last column we also give the 
corresponding effective cross- sect ions for these channels in units of fb. Of course the final 
cross-sections can be obtained only after putting in the various cuts and efficiency factors. 
These efficiency factors will have to be folded with the cross-sections given in Table [5]to get 
the final effective cross-sections for the various channels. Few clarifications on our notation 
is in order. Light charged leptons could be released in the final state through two ways: 
(i) from the decay of the heavy fermions — > l^h^ and —>■ l^H^, (ii) from the decays 
of — ^ Wi- The charged leptons released from the decays are different from those 
from in two respects. Firstly, the former carry the information on the flavor structure 
of the model as discussed in the previous sections, while the latter do not. Secondly, since 
they come from decays of the heavier E^/°, they are expected to be harder than the ones 
from decays. We refer to the charged leptons from the E^/° decays as / and the ones 
from decays as I'. The notation OSD stands for opposite sign dileptons from S^/*^ 
decays, while OSD' stands for opposite sign dileptons from decays. When we have 
one charged lepton from S^/*^ decay and an opposite sign charged lepton from decay, 
then it is denoted as 0SD{1 + I') and so on. □ 

While we provide an exhaustive list of channels for the S+S^ decay mode in Table 
ISl obviously not all of them can be effectively used at the LHC to provide smoking gun 
evidence for our two-Higgs doublet Type III seesaw model. We will highlight below a few 
of these channels which appear to be particularly promising. 

• As discussed in details before, one of the main decay channels of is ^ . 
The with mass of 40 GeV, then decays subsequently to hb pairs giving rise to a 
final state signal of a pair of opposite sign dileptons (OSD) + 4 6-jets. 

s+s- ^ trh^h^ -> trbbhl ^ 46 + OSD. 

We have seen from Table [2] that the branching ratio for is also com- 

parable. This will also produce the same collider signature of 46 + OSD. The only 
observable difference will be that the 6-jets produced from the decay will be harder 

^ We should also mention at this point that for some cases whether the charged lepton in the final state 
is a I or V can be said from the detector topology only after proper cuts have been imposed on the lepton 
transverse momentum. This will however require detailed simulation, which is outside the scope of this 
paper and will be done in an independent work. 
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SI no 


Channels 


Effective cross-section (in fb) 


1 


E+E- ^ l+l-h^h^ ^ 46 + O^D 


35.84 


2 


E+E- ^ l+l-hPhP ^2b + OSD + 2t 


3.67 


3 


E+E- ^ l+l-h^h^ OSD + 4t 


0.37 


4 


E+E- ^ l+hPH-v ^ 46 + / + 2j+ 


26.88 


5 


S+S- ^ l+h^H-v ^ 46 + 05L'(/ + /')+ i>T 


8.92 


6 


E+E" ^ l+h^R-u ^ 46 + / + r+ j6t 


4.48 


7 


E+E" ^ l+h^H'u ^ 26 + / + 2t + 2j+ j^r 


2.69 


8 


E+E- ^ l+h^H-v ^ 26 + Z + 3t+ /.t 


0.45 


9 


E+E- ^ l+hPH-v ^ 26 + + r) + 2r+ /.t 


0.9 


10 


E+E- ^ l+hPH-v ^ Z + 4t + 2j+ 


0.28 


11 


E+E- ^ l+h^H-v OSDil + /') + 4t+ 


0.04 


12 


S+S- ^ l+h^H-v ^ / + 5r+ 


0.02 


13 


E+E- ^ H+uH-u ^ 46 + 4j+ 


15.68 


14 


E+E- H+vH-v ^ 46 + 2j + i>T 


10.52 


15 


S+S" ^ H^ijR-iJ U) + 2j + -+ J)T 


5.26 


16 


E+E- ^ H+uH-jy ^ 46 + 0,SL>'+ j^t 


0.86 


17 


E+E- ^ H+uH-u ^ 46 + 2r+ 


0.43 


18 


S+S- H+uR-p ^ 46 + Ir + 1/'+ 


0.53 


19 


S+S" H+uH-u ^ 26 + 2r + 4j+ />t 


3.25 


20 


E+E- ^ H+uH-u ^ 26 + 2t + 2j + r+ /it 


2.12 


21 


E+E- H+uH-u ^ 26 + 3r + 2j+ 


1.06 


22 


S+S- ^ H^uH-u -lb + 2- + 


0.32 


23 


E+E- ^ H+uH-iy ^ 26 + 4t+ i>T 


0.08 


24 


E+E- ^ H+vH-v ^ 26 + 3t + i>T 


0.02 


25 


S+S- ^ H+uH-u ^ 4r + 4j+ 


0.15 


26 


E+E- ^ H+vH-v ^ 4r + 2j + /'+ 


0.10 


27 


E+E- H+uH-v ^ 5t + 2j+ 


0.05 


28 


E+E- ^ H+vH-v ^ 5t + r+ i)T 


0.006 


29 


E+E- ^ H+vH-v ^ 4r + i>r 


0.02 


30 


E+E- ^ H+uH'u 6r^ 


0.005 



Table 5: Effective cross-sections (in fb) for different E+E decay diannels for = 300 



as is much more massive than h^. Here and everywhere else in this section, we 
will ignore the information on the hardness of the 6-jets and present the sum of the 
cross-sections with and in the intermediate state. We should also stress that 
while we write only explicitly in the intermediate channels in the Tables, the cross- 
sections given in the final column always also include A^ as well as h^. One finds 
that the effective cross-section for this channel is 35.84 fb, which is rather high. The 
OSD released are expected to be hard, as they come from the decay of the massive 
fermions. 

Instead of decaying into bb pair, the h^s could decay into rf . If one of the decays 
into bb and the other into rr, we will get 

S+E" ^ l+rh^h^ l+rbbrf ^ 2b + OSD + 2r. 

This has an effective cross-section of 3.67 fb, which will reduce further due to the 
lower r detection efficiency. A third possibility exists where both the decay into 
rr pairs. The effective cross-section for this channel is small as can be seen from the 
Table [5l and will get smaller once the r detection efficiencies are folded. 

• The other dominant decay channel for decay is vH^. The neutrino will 
give missing energy while will decay into — >■ W^h^. The could decay 
hadronically giving 2 jets or leptonically giving either a r-jet -|- missing energy or e/fi 
lepton + missing energy. Since the lepton released in the — is important 
both for understanding the flavor structure of the mixing matrix as well as for tagging 
the channel in order to reduce the background, we consider first the case where one 
of heavy charged fermion decays into a hard charged lepton and and the other 
into a neutrino and H^. The most interesting channels in this case turn out to be: 

S+S- ^ th^R-y l+h^h^W-v 46 + / + 2j+ ^t, 

S+S" ^ l+h^H-u th^h^W'v ^Ab + l + T+ j6t, 

where for the former, the two (one from the S"*" decay and another from H~ decay) 
produce 4 6-jets, and the W~ decays produce two hadronic jets. In the latter channel, 
the W~ decays into rz/^, producing a r-jet. The effective cross-section for the former 
channel is 26.88 fb, while that for the latter is 4.48 fb. The effective cross-sections 
for the other channels with l^h^h^W~v in the intermediate states are given in Table 
[5l However, their cross-sections are smaller. 

• Finally, both the charged heavy fermions could decay through the H'^v mode. In 
this case we have the following leading order possibilities: 

S+S" ^ H^vH-y h^h^W^W'vv ^ Ab + Aj+ ^t, 
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The mode S+S^ ^ 46 + OSD'+ /px, appearing at serial number 16 in Table [5] 
could have been easy to tag as it contains 46-jets and pair of opposite sign dileptons 
coming from W"^ decay, and missing energy. However, the effective cross-section 
for this channel is relatively low. Note that none of the channels with H^vH~v in 
their intermediate state have / in their final state. For these channels therefore, it is 
impossible to say anything about the flavor structure of the model. 

8.2 S±E0 decay 

We give in Tables [6] and [3, all possible decay channels, final state configurations of parti- 
cles, and their corresponding effective cross-sections for the E^S° production and decays. 
We reiterate that the final effective cross-section after cuts will be obtained once these 
cross-sections are folded with the efficiency functions. For the leptons we follow the same 
convention for our notation as done for the previous section. 

• We begin by looking at the E^S*^ decays where l^h^ and TP —* vh^ . This 
would lead to the following final state configuration 

S^E° ^ l^h^h^v ^ 46 + /+ i>T, 

with a very large effective cross-section of 96.3 fb. This channel should be easy to 
tag. The 4 hard 6-jets come from the displaced vertices, and the lepton released is 
hard. This lepton will also carry information on the /i-r symmetric flavor structure 
of the model. Another unambiguous channel with significant effective cross-section 
coming from the l^h^h^v intermediate state is 

^ l^h%% ^ 26 + / + 2r+ i>T, 

where one of the decays into rf. 

• The other intermediate state which has very large effective cross-sections is S^S*' — > 
pH^ fhP . The would decay into W^hP , and into a lepton /' finally giving 

^ H^vh^v ^ 46 + /'+ 

with an effective cross-section of 107.4 fb. Alternatively, the W~ could instead decay 
into tVt- giving 

^ H^uh% 46 + r+ ^t, 
with effective cross-section of 53.7 fb, or decay into qq' giving 

S±SO ^ H^vh^v ^ 46 + 2j+ 
with an effective cross-section of 35.98 fb. 
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SI no 


Channels 


Eiiective cross-section (m id) 


1 


S±so ^ 46 + /+ 


96.3 


2 


E±EO ^ /±/iO/iV ^ 26 + / + 2t+ i)T 


19.7 


3 


E±E" ^ ^ / + 2r+ 


0.99 


4 


E±EO ^ ^ 46 + ir+ 


107.4 


5 


E±E° ^ H^vhPv ^ 46 + T+ i>T 


53.7 


6 


E±EO ^ H^uh^u ^ 46 + 2j+ 


35.98 


7 


E±EO ^ H^uh^u ^ 26 + 2r + 2j+ 


7.36 


8 


E±eo ^ H^uh^u ^ 26 + 2t + /'+ j^r 


2.42 


9 


E±eo ^ ^ 26 + 3r+ 


1.21 


10 


E±E" ^ H^uli^u ^ 1- + 2 j+ 


0.38 


11 


E±eo ^ H^uhPu + 4t+ 


0.12 


12 


E±eo ^ ^ 5t+ 


0.06 


13 


E±eo l^H^l^h^ ^ 46 + 2/ + 2j 


36.12 


14 


E±eo ^ l^H^l^h^ ^ 46 + 3Z(2Z + j^t 


12.04 


15 


E±eo ^ l^H^l^h^ ^Ab + 2l + 1t+ 


6.02 


1(3 


E±E" ^ I^H^I^li'^ ^ 26 + 2/ + 2- + 2 j 


7.1 


17 


E±EO ^ l^H^l^h^ ^ 26 + 3Z(2Z + Z') + 2t+ /.t 


2.4 


18 


E±eo ^ l^H^l^h^ ^ 26 + 2/ + 3r+ i6r 


1.20 


19 


E±E° ^ l^H^l^h^ ^ 2/ + 4r + 2j 


0.36 


20 


E±EO ^ l^H^l^h^ 3/(2/ + /') + 4r+ 


0.12 


21 


E±eo ^ l^H^l^h^ -^21 + 5t+ /.t 


0.06 



Table 6: Effective cross-sections (in fb) of different E^E" decay channels for Ms^ = 300 
GeV. 
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SI no 


Channels 


Effective cross-section (in fb) 


1 

± 


Z_i Z_i 7 11 1/ 11 I 7 '-tlj -\- i -\- '-tJ-\- 


±o.ou 


9 


y^y'^ — s- H^i/H^]^ Ah 4- ] 4- T -\- 9'?-l- 




Q 
O 


y^y^ H^i/H^i^ ^ Ah 4- n^Bd 4- 1'\ -i- 9i-i- 




A 


y^y'^ H^i/H^i^ ^ zL/i 4- 7' '^nd -i- /M -i- 9i-i- ^rr. 


9 1 Q 


o 


y^yo ^ n^uH^]^ ^ Ah -I- (I'^nd -I- 1'\ -I- T-\- 

Zj Zj 7 lyll I 7 ^ \J O U \l ^ I j ^ Tn^ J-'T 


1 no 


L) 


z_i z_i ^ J7 //J7 fc 7 -f- IjOUKI ^ I ] ^ / ^kT 


u.o / 


/ 


v± vo ^ ff ±7/ ff ±/T ^ 9/. HQ nfi ^ ]'] J- Ot -4- ;^ A ,, 

Zj Zj ^ O /^.O / ^ Z(/ ^ ^ / ) ^ Z/ ^ —J^ A -l 




Q 
O 


y^y^ H^uH^l^ 9h-^ T.'^lld -1- V\ -1- 9-r -1- 9i-l- 

z_i z_i ^ J J ixji t 7 -\- ijkj ly \i i J ^/ JyT 




Q 


Zj Zj ^ 11 1711 I 7 LyOi-/! t t ) OTn^ /^KT 


n 9"^ 


1 


y^y^ H^uH^i^ ^ 9h 4- r ^nd 4- ]'\ 4- '^t4- 

Zj Zj 7 £1 ty 11 I 7 Z,{J ^ 1jOLJ\1 ^ I J / 


u.uo 


1 1 

1 1 


V^VO H^uH^l^ (I'^Tld 4- l'\ 4- At 4- 944- ihrr, 
Zj Zj 7 £1 ly 11 I 7 ly yi ^ i J ^ ^1 ^ 1 JPT 


u.uu 


1 9 


y^y^ n^7/H^i^ T.'^Dd 4- 4- At 4- 9i4- hrr. 

Zj Zj 7 11 U 11 t 7 Ijtj ly [i -p i J -p ^1 ZjJ^ JPT 


n 09 




y^y^ H^uH^l^ ^ 9 ft -1- / -1- 9-r -1- Ai4- hrr, 

Zj Zj 7 11 1/ 11 I 7 ZjU ^ L ^ 1 ^ JPT 


9 78 


1 1 


v±vn ^ M^wM^l^ -^14- \t 4- \ i4- /)■,■ 


n 1 1 

V.7 . JL 1 


15 


y±^o _^ H^yHM'f ^ 46 + + 2/')+ i>T 


1.68 


16 


S±EO ^ H^vRn^^ ^2b + 3l{l + 21') + 2r+ 


0.32 


15 


S±so H^uH^P 3/(/ + 2/') + 4r+ 


0.02 


16 


S±so ^ H^uH^P ^ 46 + / + 2r+ 


0.42 


17 


E±E0 ^ H^uH^P ^ 26 + / + 4t+ 


0.08 


18 


E±eo ^ H^uH^P ^1 + 5t + 2j+ i)T 


0.04 


19 


E±EO ^ H^uH^P ^l + 6r+ 


0.004 


20 


E±E° ^ H^uH^P ^l + l' + 5r+ 


0.008 



Table 7: Effective cross-sections (in fb) of different E='=E° decay channels for = 300 
GeV. 
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• Large effective cross-section in the S^S*^ channel is also expected from the following 
decay chain 

^ l^h^l^H^ ^ 46 + 2/ + 2j, 

with effective cross-section of 36.12 fb. Both the leptons in this channel come from 
the heavy fermion decay vertices and carry the flavor information of the model. 

• S^S*^ could also decay through the intermediate states H^uH^P . This leads to 
20 possible final state particles and collider signatures. These are listed in Table [71 
However, the only one which has sizable effective cross-section is 

^ l^h^l^H^ ^Ab + l + 4j+ i>T. 

However, this channel has 4 light quark jets, which is always prone to problems with 
backgrounds. 

8.3 Backgrounds 



SI no 


Channels 


Effective cross-section 






in fb 


1 


4b +OSD 


35.84 


2 


46 + /+ 


96.3 


3 


46 + /'+ 


107.4 


4 


46 + r+ 


53.7 


5 


46 + / + 2j+ 


26.88 


6 


46 + 2/ + 2j 


36.12 


7 


46 + 3/(2/ + /')+ J^T 


12.04 



Table 8: Effective cross-sections in fb for = 300 GeV, for the most important channels 
for our model. 

In Tables [5], [6] and [7] we provided a comprehensive list of collider signature channels 
for the heavy fermions, and their corresponding effective cross-sections. In the previous 
subsection we had also discussed some of the most important channels with large effective 
cross-sections. In Table |8] we give a subset of those highlighted in sections 18.11 and 18. 2[ 
These are expected to be the most unambiguous channels, with smallest backgrounds and 
the largest signal cross-sections. In almost all channels listed in Table [HI the final collider 
signature contains 4 6-jets and a hard lepton coming from the primary heavy fermion 
decay vertex. In addition, the 4 6-jets come from the decay vertex which is significantly 
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displaced with respect to the heavy fermion decay vertex. The main source of standard 
model background for the channels with 4 6-jets and a lepton are the ttbb modes, which 
can give multiple 6-jets, leptons and missing energy. However, as mentioned many times 
before, the 6-jets come from displaced vertex and should not have any standard model 
background. Having the hard lepton in the final state further cuts down the background. 
Therefore, each of these collider channels are expected to have very little to no backgrounds. 
For a detailed signal to background analysis one requires a detailed simulation for the final 
state topology, which is outside the scope of this work. Nevertheless we add a few lines 
discussing qualitatively the possibility of backgrounds for some of the listed channels in 
Table El 

• 4b +OSD: Here the two opposite sign dileptons come from the S+S" decays. Since 
the are heavy with Ms± = 300 GeV, the leptons will be very hard and we can 
put a cut of pt ^ 100 GeV. The displaced vertex should remove all backgrounds. 

• 46 + /+ /pt: Here tibb does not directly give any background, unless one of the 
leptons from the final state is missed. However, the pt cut on the hard lepton and 
the displaced vertices should effectively remove any residual background. 

• 46 + /'+ /f)T: Here the pt cut on the lepton cannot be imposed as the lepton here 
comes from decay. However, the 4 6-jets still come from the displaced vertices 
and that should anyway take care of killing all backgrounds to a large extent. 

• 4b + l + 2jet+ pT'- The main background could again come from standard model tthb 
channels. This can also be removed by the displaced vertex and a cut of px ^ 100 
GeV for the lepton. 

• 46 + 2/ + 2j: Similar to the first case, but with 2 extra jets. 

• 46 + 3/(2/ + /')+ ^T'- Out of the 3 leptons in this channel, two are hard and one is 
relatively soft. In addition we have the hP displaced vertex. Therefore, this channel 
is expected to be absolutely background free. 

9 Conclusions 

The seesaw mechanism has remained the most elegant scheme to explain the smallness 
of the neutrino masses without having to unnaturally fine tune the Yukawa couplings to 
very small values. In the so-called Type HI seesaw, three self-conjugate {Y = 0) SU(2) 
triplet fermions are added to the standard model particle content. These exotic fermions 
are color singlets and belong to the adjoint representation of SU(2). These exotic fermions 
have Yukawa couplings with the standard model lepton doublet and the Higgs doublet. 
They also have a Major ana mass term. Once these heavy leptons are integrated out from 
the theory, we are left with a Majorana mass term for the neutrino given by the famous 
seesaw formula, where the smallness of the neutrino mass is explained by the largeness of 
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the heavy fermion mass, without having to fine tune the Yukawa couphngs to very small 
values. To generate neutrino masses rriiy ~ 0.1 eV, one requires that the heavy fermion 
mass should be ~ 10^'^ GeV. Being in the adjoint representation of SU(2), one of the most 
interesting feature of these exotic fermions is that they have gauge couplings, and therefore 
can be produced at collider experiments. The only constraint for the production of these 
particles at LHC is that their mass should be in a few 100 GeV range. However, in order to 
produce neutrino masses m,^ ~ 0.1 eV, one would then have to tune the Yukawa couplings 
to be ~ 10^^, which ruins completely the very spirit and motivation for seesaw. 

In order to circumvent this problem and preserve the motivation of the seesaw mech- 
anism, we propose an extended Type III seesaw model with two SU(2) Higgs doublets 
along with the three self-conjugate SU(2) fermion triplets. We impose an additional Z2 
symmetry such that one of the Higgs doublets, called $1, has positive charge while the 
other, called $2) has negative charge under this symmetry. In addition, we demand that 
all standard model particles have positive charge with respect to Z2 while the three new 
exotic fermion triplets are negatively charged. Therefore, $1 behaves hke the standard 
model Higgs, while $2 is coupled only to the exotic fermion triplets. As a result, the neu- 
trino mass term coming from the seesaw formula depends on the VEV of $2 {v'), while 
all other fermion masses are dependent on the VEV of $1 (v). We can therefore choose a 
value for v' such that m„ ~ 0.1 eV for exotic fermion masses ~ 100 GeV, without having 
to fine tune the Yukawa couplings to very small values. 

Another typical feature about neutrinos concern their peculiar mixing pattern which 
should be explained by the underlying theory. The current neutrino oscillation data sug- 
gests an inherent /i-r symmetry in the low energy neutrino mass matrix. It is therefore 
expected that this h-t symmetry should also exist at the high scale, either on its own or 
as a sub-group of a bigger flavor group. We imposed an exact /i-t symmetry on both the 
Yukawa coupling of the triplet fermions as well as on their Major ana mass matrix M. 
Therefore the low energy neutrino matrix m obtained after the seesaw had an in-built /x-r 
symmetry. As a result our model predicts 623 = 7r/4 and ^13 = 0. The mixing angle 612 as 
well as the mass squared differences Am^^ and Am^i are given in terms of the entries of 
Is and M. We showed how the oscillation parameters depend on the model parameters 
in Y-s and M. 

A very important and new aspect which emerged from our study concerns the mixing 
in the heavy fermion sector. It had been assumed in all past studies that the mixing 
matrices U-^, Uj^ and Uj^ which diagonalize the heavy fermion mass matrices Ms, M\jMh 
and MhM\j respectively, are almost unit matrices. However, we showed that for the case 
where M was /x-r symmetric, these matrices were highly non-trivial, and in particular had 
the last column as (0, 1/a/2, 1/\/2). We showed that this has observable consequences for 
the heavy fermion decays at the collider. We showed that flavor structure of our model 
was reflected in the pattern of heavy fermion decays into light charged leptons. These 
showed a /x-r symmetry. The state S3 ''^ decayed equally into muons and taus and almost 
never decayed into electrons. We have checked that this feature exists not only for our 
model, but for any model with an underlying flavor symmetry group that imposes yU-r 
symmetry on the heavy Majorana mass matrix M. In fact, we have made explicit checks 
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on the seesaw model proposed by Altarelli and Feruglio [36] , where ^4 was imposed as the 
flavor symmetry group. Though the proposed A4 model by Altarelli and Feruglio was a 
Type I seesaw model, it can be easily adapted to the Type III seesaw case. We found that 
the mixing matrices f/s, Ut even for that case had (0, 1/^2, 1/^2) as their last 

column. 

Having established the flavor structure of our model, we next turned to the production 
and detection of heavy fermions at LHC We discussed quantitatively and in details the 
cross-section for the heavy fermion production at LHC and their decay rates. While the 
production cross-sections for our model turned out to be same as that in all earlier calcula- 
tions done in the context of the one Higgs doublet model, the decay pattern for the heavy 
fermions in our case was found to be extremely different and unique. The /x-r permutation 
symmetry showed up in the flavor pattern of the heavy fermion decays due to the typical 
last column (0, 1/a/2, 1/"\/2) of the matrices U^,-, and t/^f. We also showed that in our 
case the decay rate of the heavy fermions was about 10^^ times larger than that found 
for the one Higgs doublet model. In fact, the heavy fermion decay rate for our model 
is 5.8 X 10^^ GeV and 4 x 10^^ GeV for 300 GeV charged and neutral heavy fermions, 
respectively. Therefore, while for the one Higgs doublet case one could attempt to look for 
displaced heavy fermion decay vertices, in our case they will decay almost instantaneously. 
We found that this tremendous decay rate came from the very fast decays of E^/*^ into 
light leptons and Higgs h'^, A'^ or H^. The very large decay rate was shown to stem from 
the very large Yukawa couplings in our model. As the Yukawa couplings are a factor of 
10^ — 10^ larger in our model, the decay rates which depend on the square of the Yukawa 
couplings are a factor lO^^-lO^^ higher. Decays into and gauge bosons in our model 
was shown to be same as in the one Higgs doublet case, and the reason explained. 

Another distinctive feature of our model appeared in the pattern of the Higgs decays. 
We showed that the smallness of the neutrino masses constrained the neutral Higgs mixing 
angle a to be very small. This resulted in a very small decay rate for the Higgs. For a 
mass of M^o = 40 GeV, the lifetime in the Higgs rest frame comes out to be about 5 
cm. This will give a displaced decay vertex in the LHC detectors, ATLAS and CMS. The 
lifetime for turned out to be small. 

Finally, we discussed in detail the expected collider signatures for our two Higgs doublet 
Type HI seesaw model with /i-r symmetry. We tabulated a comprehensive and exhaustive 
list of all possible collider signature channels for the heavy fermions at LHC. We gave the 
effective cross-sections for each of these channels. The effective cross-section for some of 
these channels were seen to be very high. We made a short-list of channels with very high 
effective cross-section and low background at LHC. This was presented in Table [H In all 
the listed channels we have 4 6-jets coming from the decay of two which have displaced 
vertices. In addition, all (except one) of them have a hard lepton in the final state coming 
from the S^/" decay. This lepton in addition to being hard, also carries information on the 
H-T symmetry of our model. These collider signature channels are hence very distinctive 
of our model and suffer from almost no standard model background. 

In conclusion, we proposed a Type HI seesaw model with large Yukawa couplings and 
triplet fermion masses light enough to be produced at the LHC. This could be achieved 
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through a unique two Higgs doublet model. The very large Yukawa couplings resulted in 
very fast decays of the heavy fermions, with a decay rate about lO^-*^ times faster than 
obtained in the earlier Type III seesaw models. We imposed a ^-r symmetry on our model 
in order to comply with the low energy neutrino oscillation data. This flavor pattern is 
reflected also in the mixing matrices of the heavy fermions, which are no longer unity, 
and which have observable consequences for the heavy fermion decays at the LHC. The 
neutrino mass constrains also the mixing angle of the neutral Higgs to be very small. This 
nearly forbids the decay of the heavy fermions into the heavier CP even neutral Higgs 
H^. Thus they decay almost always into (and A^) and More importantly, the 

small neutral Higgs mixing angle increases the hfetime of the hP. These are expected to 
live for more than 10 cm at LHC before decaying predominantly into bb, producing &-jets. 
This would be seen as a displaced decay vertex in the detector. We identified collider 
signature channels at LHC which have very large effective cross-section and almost no 
standard model background. These could be used to provide smoking gun evidence for our 
model. 
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Appendix 

A The Scalar Potential and Higgs Spectrum 

Our model has two SU(2) complex Higgs doublets $i and $2, with hypercharge Y — 1. 
The scalar potential can then be written as 

V = Ai($l$i-v')' + A2($5$2-v")' + A3(($I$i-^;') + ($^$2-i'"))' 
+A4 (($l$i)($^$2) - ($1*2) ($^$1)) + As (Re($l$2) - vv'cos^y 
+Xg (lm($l$2) - Wsin^)' , (Al) 

where 

(*i) = (°)> (*2) = (^°.,), and tan/3 = (A2) 

Recall that under the imposed Z2 symmetry, $1 carries charge +1, while $2 has —1 charge. 
Therefore, the A5 term is zero when the symmetry is exact. We will discuss shortly the 
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phenomenological consequences of this and argue in favor of a mild breaking of this Z2 
symmetry. With the scalar potential Eq. (lAip it is straightforward to obtain the Higgs mass 
matrix and obtain the corresponding mass spectrum. The physical degrees of freedoms 
contain the charged Higgs if^ and the neutral Higgs /i'', and . While and Yl^ are 
CP even, is CP odd. If we work in a simplified scenario where ^ is taken as zero, then 
it is is quite straightforward to derive the mass of the charged Higgs and the CP-odd 
Higgs . The masses are given as 

M|± = K{v^ + t;'^), and = Xe{v^ + v'^), (A3) 

respectively. The mass matrix for the neutral CP-even Higgs is 

V (4A3 + A5)W Av'\X2 + X3) + v^Xj ^ ' 

The mixing angle, obtained from diagonalizing the above matrix is given by 

tan 2a = ^ — , (A5) 

Mn - M22 ^ ' 

and the corresponding masses are 

1 



Mfio^h^ = -{Mn + M22 ± V(^ii - + 4Mf2. (A6) 

The physical Higgs are given in terms of components of $1 and $2 as follows. The neutral 
Higgs are given as 

= y2((Re<l>?-i;)cosa + (Re<l>°-i;')sma), (AT) 

/i° = v^(-(Re<l>°-t;)sina + (Re$^-t;')cosa), (A8) 

A° = V2(-Im$?sin/? + Im<l>°cos/5), (A9) 

while the charged Higgs are 

i7± = -$±sin/5 + $^cos/?. (AlO) 

The Goldstones turn out to be 

G± = $±cos/5 + $^sin/5 (All) 

G° = V2(Im$°cos/3 + Im$°sin/3). (A12) 

Recall that the requirement from small neutrino masses mj, ~ 0.1 eV constrains v' ~ 10~^ 
GeV. Therefore, for our model we get from Eqs. (lA2p and (lASp 



tan/3 ~10~^ and tan 2a ~ tan /3 ~ 10"*^. (A13) 
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One can estimate from Eq. (lA6p . that in the hmit v' <^ f , 

M^o ~ (Ai + A3)^;^ and M^o ~ X^v'^. 



(A14) 



We should point out here that in the hmit of exact Z2 symmetry, A5 = exactly, and in 
that case M^o oc ■ Since v' ~ 10~^ GeV, this would give a very tiny mass for 

the neutral Higgs h^. To prevent that, we introduce a mild explicit breaking of the Z2 
symmetry, by taking A5 7^ in the scalar potential. This not only alleviates the problem of 
an extremely light Higgs boson, it also circumvents spontaneous breaking of Z2, when the 
Higgs develop vacuum expectation value. This saves the model from complications such 
as creation of domain walls, due to the spontaneous breaking of a discrete symmetry. The 
extent of breaking of Z2 is determined by the strength of A5. Since we wish to impose only 
a mild breaking, we take A5 ~ 0.05. This gives us a light neutral Higgs mass of M° ~ 40 
GeV from Eq. ( ]A14I) . Since all other Aj ~ 1, the mass of the other CP even neutral Higgs, 
the CP odd neutral Higgs and the charged Higgs are all seen to be ~ f GeV from Eqs. 
f lA3l) and flAlil) . We will work with M% = 150 GeV, = 140 GeV and M| = 170 GeV. 

Also required are the couplings of our Higgs with the gauge bosons. This is needed in 
order to understand the Higgs decay and the subsequent collider signatures of our model. 
These are standard expressions and are well documented (see for instance [23]). One can 
check that certain couplings depend on sin a and sm{(3 — a). From Eq. ( 1A13I) we can see 
that these couplings are almost zero. Others depend on cos a and cos(/5 — a) and therefore 
large. We refer the reader to [23] for a detailed discussion on the general form for the 
coulings. 



B Appendix B: The Interaction Lagrangian 



B.l Lepton-Higgs Coupling 

The lepton Yukawa part of the Lagrangian for our two Higgs doublet model was given in 
Eq. as. 



Y 



YiJ^^V^ + + h.c] + -M,, Tr 



1 



(Bl) 



From this one can extract the individual Yukawa coupling vertex factors between two 
fermions and a Higgs. We have three generations of heavy and light neutral leptons and 
three generations of heavy and light charged leptons. In addition, we have three neutral 
and a pair of charged Higgs. The Yukawa interaction between any pair of fermions and 
a corresponding physical Higgs field can be extracted from Eq. fIBll) . We list below all 
Yukawa possible interactions in the mass basis of the particles. The vertex factors are 
denoted as Cp}^^^, where F and I are the initial and final state fermions respectively, X is 
the physical Higgs involved and L/R are for either the vertex with Pl or Pr respectively, 
where Pl and Pr are the left and right chiral projection operators respectively. Note that 
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we have suppressed the generation indices for clarity of the expressions. But the generation 
indices are imphcitly there and the vertex factors are all 3 x 3 matrices. 



= //°{I((7f '^P^ + + {I(C,5^P^ + (B2) 

+E-(cs_p,+cr#-p^)s-} 



- C>1%. = h'{lK '^P^ + C^, '^Pn)l + {/(Cf^-^L + Cf^iV«)S- + h.c} (B3) 
+E-(<4_P^ + <4_P«)E-} 



-Ctl = A°{I(Cr'^P, + + mf-"^^ (B4) 



= G°{I(Cf '^Pz. + '^P^)/ + • + Cg^^^ (B5) 



+I]'^(Cj^o^oPl + ^ j:o-£o Pr)^^} 



-<^o = /.°Mc,^:'^p, + c,^;'^p^)z.+{i7(c;jo^p^+c;;^^^^^^ (B7) 



-<^o = ^°mc5'^p, + c,i:'^p^)z. + {i7(c,^iv, + c^^^^^ (B8) 

+E°(C^o^Pl + C'^o^Pr)^'^} 

-£^0 = G°{I7(C,«:'^P^ + C£'^P«)^ + {I7(C^^^^^^ (B9) 
+EO(C^os^Pl + C'^os^Pr)^''} 
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II 

hO,L 
IT,- 
H'\L 

s-s- 

'-^11 

A°,L 
IT,- 

G0,L 
IT,- 

s-s- 



a 

IT,- 

s-s- 

U 

hO,R 
IT.- 
h'\R 

s-s- 

A0,R 
II 

A0,R 
«E- 

riA",R 

r^G°,R 
'^11 

^G°,R 
'IT.- 
,G°,R 

's-s- 



12 



12 



-^(T/iYJS'ii 

72(^11^^: 
;^(T^2>15: 
tII^ii'^'S'ii 

72(^11^^ 

^(TiiYiSi2 
^(Ti2YiSi2 



11 



12 



12 



cos a + 
cos a + 
cos a + 
sin a — 
sin a — 
sin a — 
sin P + 
sin /5 + 
sin P + 
cos (3 — 
cos /3 — 
cos P — 



72i>s5'ii sin a 



-21 

72ils5'i2 sin a 
722^s5'i2sina 
^li^S'S'ii cos a 
T^iY^Si2Cosa 

7'22^s5'l2COSa 

T'ji^S'S'ii cos/5 
T2iYj:Si2 cos/5 

722'^s5l2 COS/5 

rjiFs^n sin/5 
T2iFs5'i2 sin/5 
^22^s5'i2 sin/5 



a 
a 

c, 
a 
a 

a 



^(5iWTnCos« + 
^(5iWri2Cosa + 

^(5Lr;ri2cosa + 

^{Sl^Y^Tii^ma- 
^(5jiY'/Ti2sina - 
7|('^i2^/^i2sina - 
^(5iWTnsin/5 + 
^((5lir;Ti2sin/5 + 
^(5l2l^tTi2sin/5 + 
^(^IiF;rncos/5- 
^(^IiF;Ti2Cos/5- 
^^[Sl2Y^T,2CosP- 



S\]YIT2i sin a) 
5liys^T22sina) 
5*12^2^22 sin a) 
S\iYyT2i cos a) 
S\iY^T22 cos a) 

S'|2>s^22 COS a) 
5Lr^T2i COS/?) 

^Iiys^T22Cos/5) 

^l2>^Sr22COS/5) 

S\,YlT2i sin/?) 

5liys^r22sin/5) 

5l2rs^T22sin/5) 



Table 9: The vertex factors for (Pr) and their corresponding exact expression in terms 
of the Yukawa couphngs and mixing matrices are given in the first (third) and second 
(forth) column respectively. The vertex factors listed here are for Yukawa interactions of 
the charged leptons with neutral Higgs. 



-^fEO,.,E- = H-{l{Cr''PL + C;t~'''PR> + KCS^^''PL + CS:?''PR)^' (BIO) 

+^icS:fPL + c^:fPR)y + tF{c^:^,Pl + c^:^p«)s°} 

+h.c 



-^£o,.,E- = i/-{i(Qf •^PL + Qr'V^)z^+i(c'g;'Vi + c'g;'«p«)s° (Bii) 

+S-(Cg:'V^ + C^:fPn)u + W{Ci-4,P, + C^:§Pn)TP} 
+h.c 

The exact vertex factors Cpj^^^ for our two Higgs doublet Type III seesaw model are 
listed in Tables El [TOl [IH 

B.2 Lepton-Gauge coupling 

The lepton-gauge couplings come from the kinetic energy terms for the S fields in the 
Lagrangian. The kinetic energy terms are given as 

-A = ^zYD^^r' + L^,''\ (B12) 
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^{UlY^U^,) 




2 l'^22-'S'^12j 


r<A°,R 


i cos 
2 


^{Ul.YlU^,) 




« sin /3 1' r rT \^ rr ^ 

2 l'^21^S'^llJ 


r'G°,R 


i sin 
2 


HulY^U*,) 




2 (,tV21^S<^12j 


r^G°,R 


i sin 
2 


^UlYlU^,) 




2 l<^22^S<^12j 


riG^,R 


isin 
2 


HUI,Y^U*,) 



Table 10: The vertex factors for Pl (Pr) and their corresponding exact expression in terms 
of the Yukawa couphngs and mixing matrices are given in the first (third) and second 
(forth) column respectively. The vertex factors listed here are for Yukawa interactions of 
the neutral leptons with neutral Higgs. 



-Tl.YiUu sin P 
-Tl,YiU^2 sin P 
{^^UlYj,S,2 - Uf,Y^^S22) cos/3 
-Tl,YiU^2 sin (3 
Ti\F;f/nCos/5 
tI,Y^U^2 cos P 
{j.^UlYj:S,2 - Ul,Yt'S22) sin/5 
Tl^YiU^2Cos(5 



H-,L 

V 



,L 
G-,L 



r^G-,L 



c: 



G-,L 



r^H-,R 

H~,R 



a 



E-EO 
G-,R 



y^G ^R 

riG-,R 



a 



G-,R 
E-EO 



ij^SlY^U2i* - 4i>E*f/ii*) COS (3 
{j-^SlY^U22* - SlYs*U^2*) cos/5 

-[/i\F/Ti2sin/5 
(72'^l2>^^^^22* - Sl,Ys*U^2*) cos/3 
{j^SlY^U2i* - SlYj,*Un*) sin/3 
{j^SlYlU22* - 4i>E*f/i2*) Sin/3 

f/|ir;Ti2Cos/3 
{i^SUY^U22* - SI^Yj:*U,2*) sin/3 



Table 11: The vertex factors for Pl (Pr) and their corresponding exact expression in terms 
of the Yukawa couplings and mixing matrices are given in the first (third) and second 
(forth) column respectively. The vertex factors listed here are for Yukawa interactions of 
the charged as well as neutral leptons with charged Higgs. 
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where the first term is for heavy triplet fermion field and the second term contains the 
corresponding contributions from all standard model fields. The T^'j^ field is defined in Eq. 
(IS]). The covariant derivative is defined as 

D, = d,-V2gl^^l _^3). (B13) 



Inserting the covariant derivative in Eq. (1B12P one obtains the following interaction terms 
between leptons and gauge fields 

where the first two terms contain the neutral current interactions between and (first 
term) and between u and S*^ (second term) respectively. The last term gives the charged 
current interaction between the leptons. The neutral current interaction Lagrangian in- 
volving / and S~ is given by 

J^fc = ir{ca'''PR + 4'''PL}lZ, + {lY{c^'^PR + c^'^PL}^-Z, + h.c} 

+S^7^{44„P« + c54_P4S- Z^, (B15) 

where 



Z,R 
Cll = 




- Cu,g{Tl{T2i)} 


Z,R 




- C^„5((^2\^22), 


Z,R 


^sl{Tl,T,2) 


- Ct„5'(^22^22), 



(B16) 



Cffi- = — (-^ + S^)(5ll5i2) - C^^(4l'522), 

4_V = ^{-\^sl){S\2S^2)-c^9{S\2S22). (B17) 
The neutral current interaction Lagrangian involving the neutral leptons is given by 

C^tc = {9Cy. + 9'sy.)\vY{{UlUn)PL}iyZ^ 

Hgcn. + g's^)^TPY{{Ul2Ui2)PL}^'Z, 

Higcrn + g's^)h7Y{{Ul,Uu)PL}^'Z^ + h.c}. (B18) 
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The charged current interaction Lagrangian is given by 

dec = 9Vr{{{UlS2i) + ^iUlSu)}PL + iUlT2^)PR}lW;l (B19) 

+9T^Y{{iulS22) + -^iuls,2)}PL + (f/|;T22)p«}s-w^; 

+gTPY{{{uLS2i) + -^iUl2Sn)}PL + {UgT2i)PR}lW;^ 
+gTPr{{iUl,S22) + -^{UI,S,2)}Pl + {UgT22)PR}^-W;l + h.c 

B.3 Quark-Higgs coupling 

Finally, we discuss the the Yukawa Lagrangian for quark sector, which is given by 

-Cq = Yu^^T^MQ'.+YnA^Wj + h.c, (B20) 

where Q' is the left-handed quark doublet and u'^ and d'j^ are the right-handed "up" and 
"down" types of quark fields. Again, primes denote the flavor bases. After the electroweak 
spontaneous symmetry breaking the up and down quark mass matrices are obtained as 

Mu = Yuv (B21) 
Md = Ynv 

Note that only $i couples to both the up and down quark fields due to the imposed Z2 
symmetry, while the Yukawa couplings of $2 to quarks is forbiddeij^. However, due to the 
mixing between Higgs fields as discussed in Appendix A, all the physical Higgs particles 
would couple to the quark fields. Here we list all the interaction vertices between quarks 
and Higgs fields, which are specific to our model. The fields represents the fields in the 
mass basis. 

= _ uM^uH^ + ^ dMddH^ B22) 

V2 ^ 

= -l=^uM^uh'' - ^'^dMddh'^ (B23) 
-£f, = i^^u^^'M^uA' -i^^d^'MdA^^ (B24) 



72 


V 


1 


sina- 


72 


V 


1 

i — = 


sin(3- 
( 



lOn 



-^This is a major difference between our model and other two Higgs doublet models where the Higgs 
which couples to the neutrinos also couples to the up type quarks, while the one which couples to the 
charged leptons couples to the down type quarks. 
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£f, = -z^^MT^M^MGO + ^^^rfT^MrfdCO (B25) 



^fd = ^u{VcKMMdPR-MuVcKMPL)d + h.c (B26) 



= ^u{VcKMMdPR-M^VcKMPL)d + h.c (B27) 
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